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GPS receiver in order to reduce the consumed power of the GPS 
receiver. 

5. Claims 6, 7, 51, 61, 62, 125, 178-195, 197-203 are directed to 
calibrating a local oscillator in a mobile GPS receiver by receiving 
a reference signal from a base station via a communication link. 

6. Claims 10, 144, 196 are directed to the position determination 

of a remote unit, whereby the measurements are partially made at the 
base station, Doppler and ephemeris measurements, and partially at 
the remote unit, pseudorange measurements; the measurements of the 
remote unit are transmitted to the base station, where position 
calculation takes place. 



INTERNATIONAL SEARCH REPORT 

Information on patent family members 


(ni ionai Application No 

PCT/US 96/16161 


Patent document 
cited in search report 


Publication 
dale 


Patent family 
member (s) 


Publication 
date 



US 4998111 A 05-03-91 NONE 



US 4445118 


A 


24-04-84 


NONE 




GB 2273218 


A 
n 








04-07-95 








DE 


4340954 A 


09-Q6-94 








JP 


6213992 A 


05-08-94 


US 4785463 


A 


15-11-88 


AU 


6597086 A 


24-03-87 








CA 


1260120 A 


26-09-89 








EP 


G235285 A 


nQ.nq.87 

\jj of 








WO 


8701540 A 


1P-G3-87 

XC UJ Of 


WO 9428434 


A 

A 


08-1^-34 


US 


5418538 A 


cj 


US 5491486 


A 


13-02-96 


CA 


2163628 A 


02-11-95 






EP 


0705441 A 


10-04-96 








JP 


8512163 T 


17-12-96 








W0 


9529410 A 


02-11-95 


EP 0545636 


A 


09-06-93 


US 


5379224 A 


03-01-95 






AU 


667492 B 


28-03-96 








AU 


2965792 A 


03-06-93 








CA 


2083343 A 


30-05-93 








FI 


925403 A 


30-05-93 








JP 


5256927 A 


08-10-93 



US 5448773 A 05-09-95 NONE 



EP 0447978 


A 


25-09-91 


JP 3269385 A 
DE 69111275 D 
DE 69111275 T 
US 5222245 A 


29-11-91 
24-08-95 
08-02-96 
22-06-93 


US 5365450 


A 


15-11-94 


AU 6226794 A 
CA 2152052 A 
W0 9415412 A 
US 5422813 A 


19-07-94 
07-07-94 
07-07-94 
06-06-95 



US 5245634 A 14-09-93 CA 2090740 A,C 24-09-93 



CORRECTED 
VERSION* 

PCT 



WORLD INTELLECTUAL PROPERTY ORGANIZATION 
International Bureau 




(51) International Patent Classification 0 : 




(11) International Publication Number: 


WO 97A4049 


G01S L/04, 5/14, 1/02 


A3 


(43) International Publication Date: 


17 April 1997 (17.04.97) 



(21) International Application Number: PCT/US96/16161 

(22) International Filing Date: 8 October 1996 (08.10.96) 



(30) Priority Data: 

60/005.318 
08/612 t 582 
08/612,669 
08/613,966 



9 October 1995 (09.10.95) 
8 March 1996 (08.03.96) 
8 March 1996 (08.03.96) 
8 March 1996 (08.03.96) 



(71) Applicant: SNAPTRACK, INC. [US/US]; Suite 250, 
Moorpark Avenue, San Jose, CA 951 17 (US). 



US 
US 
US 

us 



4040 



(72) Inventor: KRASNER, Norman, F. 
Carlos, CA 94070 (US). 



117 Coventry Court, San 



(74) Agents: SCHELLER, James. C. et al.; Blakeiy, Sokoloff, 
Taylor & Zafman, 7th floor. 12400 Wilshire Boulevard, Los 
Angeles, CA 90025 (US). 



(81) Designated States: AL, AM. AT, AT (Utility model), AU..AZ, 
BA, BB. BG, BR. BY. CA, CH, CN, CU, CZ, CZ (Utility 
model), DE, DE (Utility model), DK, DK (Utility model)," 
EE, EE (Utility model), ES, FI, FI (Utility model). GB, GE, 
HU IL, IS, JP. KE, KG, KP, KR, KZ, LC, LK. LR. LS, 
LT LU, LV, MD, MG, MK. MN, MW, MX. NO, NZ, PL. 
PT, RO, RU, SD. SE, SG, SI, SK, SK (Utility model), TJ, 
TM, TR, TT. UA, UG. UZ. VN. ARIPO patent (KE, LS. 
MW, SD, SZ, UG), Eurasian patent (AM, AZ, BY, KG. KZ. 
MD RU. TJ. TM). European patent (AT. BE. CH. DE, DK, 
ES FI FR. GB. GR, IE, IT. LU, MC NL, PT. SE). OAPI 
patent (BF. BJ, CF, CG, CI, CM, GA, GN, ML, MR, NE. 
SN, TD, TG). 



Published 

With international search report. 

Before the expiration of the time limit for amending the 
claims and to be republished in the event of the receipt of 
amendments. 

(*&) Date of publication of the international search report: 

1 19 June 1997 (19.06.97) 



(54) Title: GPS RECEIVER AND METHOD FOR PROCESSING GPS SIGNALS 



Mobile or Remote Unit 



GPS Antenna 



26 



(57) Abstract 
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downconverter coupled to 
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received GPS signals to 
an intermediate frequency 
(IF); a digitizer ..coupled 
to the downconverter and 
sampling the IF GPS 
signals at a predetermined 
rate to produce sampled 
IF GPS signals; a memory 
coupled to the digitizer 
storing the sampled IF GPS 
signals (a snapshot of GPS 
signals); and a digital signal 
processor (DSP) coupled to 
the memory and operating 
under stored instructions 
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features and includes, in another embodiment the capability to correct for errors in its local oscillator which is used to sample the GPS 
signals. The calculation speed of pseudoranges, and sensitivity of operation, is enhanced by the transmission of the Doxwler freouencv 
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GPS RECEIVER AND METHOD FOR PROCESSING 
GPS SIGNALS 

BACKGROUND OF THF INWNTp^ 

RELATED APPLICATIONS 

This application is related to two patent applications filed by the same 
inventor on the same date as this application; these two applications are: An 
Improved GPS Receiver Utilizing a Communication Link (Serial No. 08/612,582, 
filed on March 8, 1996); An Improved GPS Receiver Having Power Management 
(Serial No. 08/613,966, filed on March 8, 1996). 

This application is also related to and hereby claims the benefit of the filing 
date of a provisional patent application by the same inventor, Norman F. Krasner, 
which application is entitled Low Power, Sensitive Pseudorange Measurement 
Apparatus and Method for Global Positioning Satellites Systems, Serial No. 
60/005,318, filed October 9, 1995. 

A portion of the disclosure of this patent document contains 
material which is subject to copyright protection. The copyright owner has 
no objection to the facsimile reproduction by anyone of the patent 
document or the patent disclosure, as it appears in the Patent and 
Trademark Office patent file or records, but otherwise reserves all 
copyright rights whatsoever. 

1. FIELD OF THF INVENTION 

The present invention relates to receivers capable of determining 
position information of satellites and, in particular, relates to such receivers 
which find application in global positioning satellite (GPS) systems. 

2. BACKGROUND ART 

GPS receivers normally determine their position by computing 
relative times of arrival of signals transmitted simultaneously from a 
multiplicity of GPS (or NAVSTAR) satellites. These satellites transmit, as 
pan of their message, both satellite positioning data as well as data on clock 
timing, so-called "ephemeris" data. The process of searching for and 
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acquiring GPS signals, reading the ephemeris data for a multiplicity of 
satellites and computing the location of the receiver from this data is time 
consuming, often requiring several minutes. In many cases, this lengthy 
processing time is unacceptable and, furthermore, greatly limits battery life 
in micro-miniaturized portable applications. 

Another limitation of current GPS receivers is that their operation is 
limited to situations in which multiple satellites are clearly in view, without 
obstructions, and where a good quality antenna is properly positioned to 
receive such signals. As such, they normally are unusable in portable, body 
mounted applications; in areas where there is significant foliage or building 
blockage; and in in-building applications. 

There are two principal functions of GPS receiving systems: (1 ) 
computation of the pseudoranges to the various GPS satellites, and (2) 
computation of the position of the receiving platform using these 
pseudoranges and satellite timing and ephemeris data. The pseudoranges are 
simply the time delays measured between the received signal from each 
satellite and a local clock. The satellite ephemeris and timing data is 
extracted from the GPS signal once it is acquired and tracked. As stated 
above, collecting this information normally takes a relatively long time (30 
seconds to several minutes) and must be accomplished with a good received 
signal level in order to achieve low error rates. 

Virtually all known GPS receivers utilize correlation methods to 
compute pseudoranges. These correlation methods are performed in real 
time, often with hardware correlators. GPS signals contain high rate 
repetitive signals called pseudorandom (PN) sequences. The codes available 
for civilian applications are called C/A codes, and have a binary phase- 
reversal rate, or "chipping" rate, of 1.023 MHz and a repetition period of 
1023 chips for a code period of 1 msec. The code sequences belong to a 
family known as Gold codes. Each GPS satellite broadcasts a signal with a 
unique Gold code. 

For a signal received from a given GPS satellite, following a 
downconversion process to baseband, a correlation receiver multiplies the 
received signal by a stored replica of the appropriate Gold code contained 
within its local memory, and then integrates, or lowpass filters, the product 
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in order to obtain an indication of the presence of the signal. This process is 
termed a "correlation" operation. By sequentially adjusting the relative 
timing of this stored replica relative to the received signal, and observing the 
correlation output, the receiver can determine the time delay between the 
received signal and a local clock. The initial determination of the presence of 
such an output is termed "acquisition." Once acquisition occurs, the process 
enters the "tracking" phase in which the timing of the local reference is 
adjusted in small amounts in order to maintain a high correlation output. 
The correlation output during the tracking phase may be viewed as the GPS 
signal with the pseudorandom code removed, or, in common terminology, 
"despread." This signal is narrow band, with bandwidth commensurate 
with a 50 bit per second binary phase shift keyed data signal which is 
superimposed on the GPS waveform. 

The correlation acquisition process is very time consuming, 
especially if received signals are weak. To improve acquisition time, most 
GPS receivers utilize a multiplicity of correlators (up to 12 typically) which 
allows a parallel search for correlation peaks. 

Some prior GPS receivers have used FFT techniques to determine 
the Doppler frequency of the received GPS signal. These receivers utilize 
conventional correlation operations to despread the GPS signal and provide 
a narrow band signal with bandwidth typically in the range of 10 kHz to 30 
kHz. The resulting narrow band signal is then Fourier analyzed using FFT 
algorithms to determine the carrier frequency. The determination of such a 
carrier simultaneously provides an indication that the local PN reference is 
adjusted to the correct phase of the received signal and provides an accurate 
measurement of carrier frequency. This frequency may then be utilized in 
the tracking operation of the receivers. 

U.S. Patent No. 5,420,592 to Johnson discusses the use of FFT 
algorithms to compute pseudoranges at a central processing location rather 
than at a mobile unit. According to that method, a snapshot of data is 
collected by a GPS receiver and then transmitted over a data link to a remote 
receiver where it undergoes FFT processing. However, the method 
disclosed therein computes only a single forward and inverse Fast Fourier 
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Transform (corresponding to four PN periods) to perform the set of 
correlations. 

As will be evident from the following description of the present 
invention, higher sensitivity and higher processing speed can be achieved by 
performing a large number of FFT operations together with special 
preprocessing and postprocessing operations. 

In this patent the terms correlation, convolution and matched filtering 
are often utilized. The term "correlation" when applied to two series of 
numbers means the term by term multiplication of corresponding members 
of the two series followed by the summation of the series. This is 
sometimes referred to as "serial correlation" and results in an output that is a 
single number. In some circumstances, a succession of correlation 
operations are performed on successive groups of data. 

The term "convolution" as applied to two series of numbers is the 
same as that commonly used in the an and is equivalent to a filtering of the 
second series of length m with a filter, corresponding to the first series, 
having an impulse response of length n. The result is a third series of length 
m+n- 1. The term "matched filtering" refers to a convolution, or filtering, 
operation in which the aforementioned filter has an impulse response which 
is the time-reversed complex conjugate of the first series. The term "fast 
convolution" is utilized to indicate a series of algorithms that computes the 
convolution operation in an efficient manner. 

Some authors utilize the terms correlation and convolution 
interchangeably; for clarity, however, in this patent, the term correlation 
always refers to the serial correlation operation described above. 

SUMMARY 

One embodiment of the present invention provides a method for 
determining the position of a remote GPS receiver by transmitting GPS satellite 
information, including Doppler, to the remote unit or mobile GPS unit from a 
basestation via a data communication link. The remote unit uses this information 
and received GPS signals from in view satellites to subsequently compute 
pseudoranges to the satellites. The computed pseudoranges are then transmitted to 
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the basestation where the position of the remote unit is calculated. Various 
embodiments of apparatuses which can perform this method are also described. 

Another embodiment of the present invention provides a GPS receiver 
having an antenna for receiving GPS signals from in view satellites; and a 
downconvener for reducing the RF frequency of the received GPS signals to an 
intermediate frequency (IF). The IF signals are digitized and stored in memory for 
later processing in the receiver. This processing typically is accomplished, in one 
embodiment of the invention, using a programmable digital signal processor which 
executes the instructions necessary to perform fast convolution (e.g. FFT) 
operations on the sampled IF GPS signals to provide pseudorange information. 
These operations also typically include preprocessing (prior to fast convolution) 
and post processing (after fast convolution) of stored versions of the GPS signals 
or processed and stored versions of the GPS signals. 

Yet another embodiment of the present invention provides a method of 
power management for a GPS receiver and also provides a GPS receiver having 
power management features. Power dissipation is reduced over prior systems by 
receiving GPS signals from in view satellites; buffering these signals; and then 
turning off the GPS receiver. Other power management features are described. 

BRIEF DESCRIPTION OF T HE DRAWINGS 

The present invention is illustrated by way of example and not limitation in 
the figures of the accompanying drawings in which references indicate similar 
elements and in which: 

Figure 1A is a block diagram of the major components of a remote or 
mobile GPS receiving system utilizing the methods of the present invention, 
and shows data links that may exist between a basestation and the remote. 

Figure IB is a block diagram of an alternative GPS mobile unit. 

Figure 1C is a block diagram of another alternative GPS mobile unit. 

Figures 2 A and 2B provide two alternatives for the RF and IF 
portions of a receiver which is an embodiment of the present invention. 

Figure 3 shows a flow chart of the major operations (e.g. software 
operations) performed by the programmable DSP processor in accordance 
with the methods of the present invention. 
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Figure 4 illustrates the signal processing waveforms at various 
stages of processing according to the methods of the present invention. 

Figure 5A illustrates a basestation system in one embodiment of the 
present invention. 

Figure 5B illustrates a basestation system in an alternative 
embodiment of the present invention. 

Figure 6 illustrates a GPS mobile unit having, according to one 
aspect of the present invention, local oscillator correction or calibration. 

Figure 7 is a flow chart which shows a power management method 
for a mobile unit according to one embodiment of the present invention. 

DETAILED DESCRIPTION OF THE INVENTION 

This invention concerns apparatuses and methods for computing the 
position of a mobile, or remote, object in a manner that results in the remote 
hardware having very low power dissipation and the ability to operate with 
very low received signal levels. That is, power consumption is reduced 
while receiver sensitivity is increased. This is made possible by the 
implementation of the remote receiving functions, as shown in Figure 1 A, as 
well as the transmission of Doppler information from a separately located 
basestation 10 to the remote or GPS mobile unit 20. 

It should be noted that pseudoranges may be used to compute the 
remote's geographical position in many different ways. Three examples 
arc: 

1 . Method 1: By re-transmitting the satellite data messages to the 
remote 20 from the basestation 10, the remote 20 may combine this 
information with the pseudorange measurements to compute its 
position. See, for example, U.S. patent No. 5,365,450, which is 
incorporated herein by reference. Typically, the remote unit 20 
performs the computation of position in the remote 20. 

2 . Method 2: The remote 20 may gather the satellite ephemeris data 
from the reception of GPS signals in the normal manner that is 
commonly practiced in the art. This data, which typically is valid 
for one to two hours, may be combined with pseudorange 
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measurements to complete, typically in the remote unit, the position 
calculation. 

3 . Method 3: The remote 20 may transmit over a communications link 
16 the pseudoranges to the basestation 10 which can combine this 
information with the satellite ephemeris data to complete the 
position calculation. See, for example, U.S. Patent No. 
5,225,842, which is incorporated herein by reference. 
In approaches (or Methods) 1 and 3, it is assumed that the 
basestation 10 and remote 20 have a common view of all satellites of 
interest and are positioned close enough to one another to resolve a time 
ambiguity associated with the repetition rate of the GPS pseudorandom 
codes. This will be met for a range between basestation 10 and remote 20 
of 1/2 rimes the speed of light times the PN repetition period ( 1 
millisecond), or about 150 km. 

In order to explain the current invention, it is assumed that method 
3 is utilized to complete the position calculation. However, upon review of 
this Specification, it will be appreciated by those skilled in the an that the 
various aspects and embodiments of the present invention could be used 
with any of the above three Methods as well as other approaches. For 
example, in a variation of Method 1, satellite data information such as data 
representative of satellite ephemeris may be transmitted by a basestation to 
a remote unit, and this satellite data information may be combined with 
pseudo ranges, computed according to the present invention from buffered 
GPS signals, to provide a latitude and longitude (and in many cases also an 
altitude) for the remote unit. It will be appreciated that the position 
information received from the remote may be limited to latitude and 
longitude or may be extensive information which includes latitude, 
longitude, altitude, velocity and bearing of the remote. Moreover, the local 
oscillator correction and/or the power management aspects of the present 
invention may be utilized in this variation of Method 1. Furthermore, 
Doppler information may be transmitted to the remote unit 20 and utilized 
by the remote unit 20 in accordance with aspects of the present invention. 

Under Method 3, the basestation 10 commands the remote 20 to 
perform a measurement via a message transmitted over a data 
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communications link 16 as shown in Figure 1 A. The basestation 10 also 
sends within this message Doppler information for the satellites in view, 
which is a form of satellite data information. This Doppler information 
typically is in the format of frequency information, and the message will 
typically also specify an identification of the particular satellites in view or 
other initialization data. This message is received by a separate modem 22 
that is pan of the remote unit 20, and it is stored in a memory 30 coupled to 
a low-power microprocessor 26. The microprocessor 26 handles data 
information transfer between the remote unit processing elements 32-48 
and the modem 22, and it controls power management functions within the 
remote receiver 20, as will be evident in the subsequent discussion. 
Normally, the microprocessor 26 sets most or all remote unit 20's 
hardware to a low power, or power down, state, except when the 
pseudorange and/or other GPS calculations are being performed, or when 
an alternative source of power is available. However, the receiver portion 
of the modem is at least periodically turned on (to full power) to determine 
if the basestation 10 has sent a command to determine the remote's 
position. 

This above-mentioned Doppler information is very short in 
duration since the required accuracy of such Doppler information is not 
high. For example, if 10 Hz accuracy were required and the maximum 
Doppler is approximately ±7kHz, then an 1 1 bit word would suffice for 
each satellite in view. If 8 satellites were in view, then 88 bits would be 
required to specify all such Dopplers. The use of this information 
eliminates the requirement for the remote 20 to search for such Doppler, 
thereby reducing its processing time by in excess of a factor of 10. The 
use of the Doppler information also allows the GPS mobile unit 20 to 
process more quickly a sample of GPS signals and this tends to reduce the 
amount of time for which the processor 32 must receive full power in order 
to compute a position information. This alone reduces the power 
consumed by the remote unit 20 and contributes to improved sensitivity. 
Additional information may also be sent to the remote 20, including the 
epochs of the data in the GPS message. 
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The received data link signal may utilize a precision carrier 
frequency. The remote receiver 20 may employ, as shown in Figure 6 
which is described below, an automatic frequency control (AFC) loop to 
lock to this carrier and thereby further calibrate its own reference oscillator. 
A message transmission time of 10 msec, with a received signal to noise 
ratio of 20 dB, will normally allow frequency measurement via an AFC to 
an accuracy of 10 Hz or better. This will typically be more than adequate 
for the requirements of the present invention. This feature will also 
enhance the accuracy of the position calculations which are performed, 
either conventionally or using the fast convolution methods of the present 
invention. 

In one embodiment of the invention, the communication link 16 is a 
commercially available narrow bandwidth radio frequency communication 
medium, such as a two-way pager system. This system may be used in 
embodiments where the amount of data transmitted between the remote 20 
and basestation 10 is relatively small. The amount of data required for the 
transmission of Doppler and other data (e.g. initialization data such as the 
identities of the satellites in view) is relatively small and similarly the 
amount of data required for the position information (e.g.. pseudoranges) 
is relatively small. Consequently, narrowband systems are adequate for 
this embodiment. This is unlike those systems which require the 
transmission of large amounts of data over a short period of time; these 
systems may require a higher bandwidth radio frequency communication 
medium. 

Once the remote 20 receives a command (e.g., from the basestation 
10) for GPS processing together with the Doppler information, the 
microprocessor 26 activates the RF to IF Convener 42, Analog to Digital 
Converter 44 and Digital Snapshot Memory 46 via a Battery and Power 
Regulator and Power Switches circuit 36 (and controlled power lines 21a, 
21b, 21c and 2 Id) thereby providing full power to these components. 
This causes the signal from the GPS satellite which is received via antenna 
40 to be downconverted to an IF frequency, where it subsequently 
undergoes digitization. A contiguous set of such data, typically 
corresponding to a duration of 100 milliseconds to 1 second (or even 
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longer), is then stored in a Snapshot Memory 46. The amount of data 
stored may be controlled by the microprocessor 26 such that more data 
may be stored in the memory 46 (to obtain better sensitivity) in those 
situations when conserving power is not as important as obtaining better 
sensitivity, and less data may be stored in those situations when 
conservation of power is more important than sensitivity. Typically, 
sensitivity is more important when the GPS signals may be obstructed 
partially, and power conservation is less important when a copious power 
supply (e.g. a car battery) is available. The addressing of this memory 46 
to store this data is controlled by a Field Programmable Gate Array 
integrated circuit 48. Downconversion of the GPS signal is accomplished 
using a frequency synthesizer 38 which provides local oscillator signal 39 
to the converter 42 as discussed further below. 

Note that all this time (while the snapshot memory 46 is being filled 
with the digitized GPS signals from the in view satellites) the DSP 
microprocessor 32 may be kept in a low power state. The RF to IF 
Convener 42 and Analog to Digital Converter 44 are typically only turned 
on for a short period of time, sufficient to collect and store the data 
required for pseudorange calculation. After the data collection is complete, 
these converter circuits are turned off or power is otherwise reduced via 
controlled power lines 21b and 21c (while the memory 46 continues to 
receive full power), thus not contributing to additional power dissipation 
during the actual pseudorange calculation. The pseudorange calculation is 
then performed using, in one embodiment, a general purpose, 
programmable digital signal processing IC 32 (DSP), as exemplified by a 
TMS320C30 integrated circuit from Texas Instruments. This DSP 32 is 
placed in an active power state by the microprocessor 26 and the circuit 36 
via controlled power line 21e prior to performing such calculations. 

This DSP 32 differs from others used in some remote GPS units in 
that it is general purpose and programmable, as compared to specialized 
custom digital signal processing ICs. Furthermore, the DSP 32 makes 
possible the use of a Fast Fourier Transform (FFT) algorithm, which 
permits very rapid computation of the pseudoranges by performing rapidly 
a large number of correlation operations between a locally generated 
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reference and the received signals. Typically, 2046 such correlations are 
required to complete the search for the epochs of each received GPS 
signal. The Fast Fourier Transform algorithm permits a simultaneous and 
parallel search of all such positions, thus speeding the required 
computation process by a factor of 10 to 100 over conventional 
approaches. 

Once the DSP 32 completes its computation of pseudoranges for 
each of the in view satellites, it transmits, in one embodiment of the 
invention, this information to the microprocessor 26 via interconnect bus 
33. At this time the microprocessor 26 may cause the DSP 32 and memory 
46 to again enter a low power state by sending an appropriate control 
signal to the Battery and Power Regulator circuit 36. Then, the 
microprocessor 26 utilizes a modem 22 to transmit the pseudorange data 
over a data link 16 to the basestation 10 for final position computation. In 
addition to the pseudorange data, a time tag may simultaneously be 
transmitted to the basestation 10 that indicates the elaspsed time from the 
initial data collection in the buffer 46 to the time of transmission of the data 
over the data link 16. This time tag improves the capability of the 
basestation to compute position calculation, since it allows the computation 
of the GPS satellite positions at the time of data collection. As an 
alternative, in accordance with Method 1 above, the DSP 32 may compute 
the position (e.g. latitude, longitude or latitude, longitude and altitude) of 
the remote unit and send this data to the microprocessor 26, which 
similarly relays this data to the basestation 10 via the modem 22. In this 
case the position computation is eased by the DSP maintaining the elapsed 
time from the reception of satellite data messages to the time at which the 
buffer data collection begins. This improves the capability of the remote 
unit to compute position calculation, since it allows the computation of the 
GPS satellite positions at the time of data collection. 

As shown in Figure 1 A, modem 22, in one embodiment, utilizes a 
separate antenna 24 to transmit and receive messages over data link 16. It 
will be appreciated that the modem 22 includes a communication receiver 
and a communication transmitter which are alternatively coupled to the 
antenna 24. Similarly, basestation 10 may use a separate antenna 14 to 
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transmit and receive data link messages, thus allowing continuous 
reception of GPS signals via GPS antenna 12 at the basestation 10. 

It is expected, in a typical example, that the position calculations in 
the DSP 32 will require less than a few seconds of time, depending upon 
the amount of data stored in the digital snapshot memory 46 and the speed 
of the DSP or several DSPs. 

It should be clear from the above discussion that the remote unit 20 
need only activate its high power consumption circuits for a small fraction 
of time, if position calculation commands from the basestation 10 are 
infrequent It is anticipated, in at least many situations, that such 
commands will result in the remote equipment being activated to its high 
power dissipation state only about 1 % of the time or less. 

This then allows battery operation for 100 times the length of time 
that would otherwise be possible. The program commands necessary for 
the performance of the power management operation are stored in 
EEPROM 28 or other suitable storage media. This power management 
strategy may be adaptable to different power availability situations. For 
example, when prime power is available the determination of position may 
occur on a continuing basis. 

As indicated above, the digital snapshot memory 46 captures a 
record corresponding to a relatively long period of time. The efficient 
processing of this large block of data using fast convolution methods 
contributes to the ability of the present invention to process signals at low 
received levels (e.g., when reception is poor due to partial blockage from 
buildings, trees, etc.). All pseudoranges for visible GPS satellites are 
computed using this same buffered data. This provides improved 
performance relative to continuous tracking GPS receivers in situations 
(such as urban blockage conditions) in which the signal amplitude is 
rapidly changing. 

A slightly different implementation exhibited in Figure IB 
dispenses with the microprocessor 26 and its peripherals (RAM 30 and 
EEPROM 28) and replaces its functionality with additional circuitry 
contained within a more complex FPGA (field programmable gate array) 
49. In this case the FPGA 49, a low power device, serves to wake-up the 
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DSP 32a chip upon sensing activity from the modem 22 through 
interconnect 19. Interconnect 19 couples the modem to the DSP 32a and to 
the FPGA 19. The DSP chip 32a, when awake, directly transmits and 
receives data from the modem. The DSP 32a also performs power control 
operations through its interconnect 18 which is coupled to the Battery and 
Power Regulator and Switches 36 to provide power on/off commands to 
the circuit 36. The DSP 32a selectively powers on or reduces power to 
different components, according to a power management method such as 
that shown in Figure 7, through the power on/off commands provided by 
interconnect 18 to the circuit 36. The circuit 36 receives these commands 
and selectively provides power (or reduces power) to the different 
components. The circuit 36 wakes up the DSP 32a via interconnect 17. 
The circuit 36 selectively provides power to the different components by 
selectively switching power through selected ones of controlled power 
lines 21a, 21b, 21c, 2 Id and 2 If. Thus, for example, to provide power to 
the converter 42 and the converter 44, power is provided through lines 21b 
and 21c to these conveners. Similarly, power to the modem is supplied 
through controlled power line 2 If. 

A low frequency crystal oscillator 47 is coupled to the memory and 
power management FPGA 49. In one embodiment, the memory and 
power management FPGA 49 contains a low power timer which includes 
the low frequency oscillator 47. When the FPGA 49 timer expires, the 
FPGA 49 sends a wake up signal to the DSP 32a through interconnect 17, 
and the DSP 32a can then wake up other circuitry by providing power 
on/off commands to the Battery and Power Regulator and Power Switches 
circuit 36. The other circuitry is powered, through the controlled power 
lines 21a, 21b, 21c, 21d and 21f under control of the circuit 36, in order to 
perform a positioning operation (e.g. determine a position information 
such as a pseudorange or a (latitude and longitude). Following the 
positioning operation, the DSP 32A resets the FPGA timer and reduces 
power to itself, and the circuit 36 also reduces power to the other 
components, in accordance with the method shown in Figure 7. It will be 
appreciated that a battery or a plurality of batteries will provide power for 
all power controlled circuits through controlled power lines which are 
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controlled by the memory and power management FPGA 49 and the DSP 
32a. It will also be appreciated that, rather than directly reducing power by 
controlling power lines (such as 21b) to a component, the power 
consumed by a component may be reduced by signaling to the component 
(as in the case of DSP 32a via interconnect 17 in Figure IB) to reduce 
power or wake up to full power, this is often possible when a component, 
such as an integrated circuit, has an input for controlling the power state of 
the component, and the component has the necessary internal logic for 
controlling power consumption (e.g. logic for reducing power to various 
logical blocks of the component). The memory and power management 
FPGA 49 provides memory control and management, including addressing 
operations when data is being stored into the memory 46 from the 
conveners 44 or when the DSP component 32a is reading data from the 
memory 46. The FPGA 49 may also be controlling other memory 
functions such as memory refresh if necessary. 

Figure 1C shows another embodiment according to the present 
invention of a GPS mobile unit which contains many of the same 
components as the GPS mobile units shown in Figures 1 A and IB. In 
addition, the GPS mobile unit shown in Figure 1C includes power 
regulators 77 which are coupled to receive power from a plurality of 
batteries 81 as well as an optional external power source input 83 and solar 
cells 79. The power regulator 77 provides power for all circuits under 
control of the controlled power lines which are managed by the DSP chip 
32a and the memory and power management FPGA 49 shown in Figure 
1C The solar cell 79 may recharge, using conventional recharging 
technology, those batteries. The solar cells 79 may also provide power to 
the GPS mobile unit in addition to recharging the batteries. In the 
embodiment shown in Figure 1C, the FPGA 49 provides a wake-up signal 
over interconnect 75 to the DSP chip 32a; this signal causes the DSP chip 
to return to full power to perform the various functions described for the 
DSP chip 32a. The DSP chip may also be activated to full power state via 
an external command from the modem 22 which is coupled directly to the 
DSP chip via interconnect 19. 
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Figure 1C also shows a feature of the present invention which 
allows the GPS mobile unit to trade off sensitivity for power conservation. 
As described herein sensitivity of the GPS mobile unit may be increased by 
increasing the amount of buffered GPS signals which are stored in the 
memory 46. This is done by acquiring and digitizing more GPS signals 
and storing this data in the memory 46. While this increased buffering 
causes more power consumption, it does improve the sensitivity of the 
GPS mobile unit. This increased sensitivity mode may be selected by a 
power mode switch 85 on the GPS unit which is coupled to bus 19 to 
provide a command to the DPS chip 32a to enter an increased sensitivity 
mode. This power mode switch 85 may alternatively be caused to send a 
command to the DSP 32a chip to conserve more power and provide less 
sensitivity by acquiring a smaller snapshot of the GPS signals and thereby 
storing a smaller amount of GPS signals in the memory 46. It will be 
appreciated that this power mode selection may also occur through a signal 
sent from the basestation to the modem 22 which then communicates this 
command via interconnect 19 to the DSP chip 32a. 

A representative example of an RF to IF frequency converter and 
digitizing system for the mobile GPS unit is shown in Figure 2A. The 
input signal at 1575.42 MHz is passed through a bandlimiting filter (BPF) 
50 and low noise amplifier (LNA) 52 and sent to a frequency conversion 
stage. The local oscillator (LO) 56 used in this stage is phase locked (via 
PLL 58) to a 2.048 MHz (or harmonic thereof) temperature compensated 
crystal oscillator (TCXO) 60. In a preferred implementation, the LO 
frequency would be 1531.392 MHz, which is 2991 x 0.512 MHz. The 
resulting IF signal is then centered at 44.028 MHz. This IF is desirable 
due to the availability of low cost components near 44 MHz. In particular, 
surface acoustic wave filters (SAW), which are utilized in abundance in 
television applications, are readily available. Of course, other bandlimiting 
devices could be used instead of SAW devices. 

The received GPS signal is mixed with the LO signal in mixer 54 to 
produce the IF signal. This IF signal is passed through a SAW filter 64, 
for precision bandlimiting to 2 MHz bandwidth, and then sent to an l/Q 
down-convener 68, which translates the signal to near baseband (4 kHz 
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center frequency nominally). The local oscillator frequency for this 
downconverter 68 is derived from the 2.048 MHz TCXO 60 as the 43rd 
harmonic of 1.024 MHz, that is 44.032 MHz. 

The I/Q downconverter 68 is generally commercially available as an 
RF component. It typically consists of two mixers and lowpass filters. In 
such instances, the input ports of one mixer are fed with the IF signal and 
the LO signal and the input pons to the other mixer are fed with the same 
IF signal and the LO signal phase shifted by 90°. The outputs of the two 
mixers are lowpass filtered to remove feedthrough and other distortion 
products. 

As shown in Figure 2A, amplifiers 62 and 66 may be used before 
and after the bandlimiting operation as required. 

The two outputs of the I/Q downconverter 68 are sent to two 
matched A/D conveners 44 which sample the signals at 2.048 MHz. An 
alternative implementation replaces the A/D conveners 44 with comparators 
(not shown), each of which outputs a two- valued (one-bit) sequence of 
data in accordance with the polarity of the incoming signal. It is well 
known that this approach results in a loss of approximately 1.96 dB in 
receiver sensitivity relative to a multilevel A/D convener. However, there 
may be substantial cost savings in use of a comparator vs. A/D converters, 
as well as in the reduced memory requirement in the following snapshot 
memory 46. 

An alternative implementation of the downconverter and A/D system is 
shown in Figure 2B which utilizes a bandpass sampling method. The TCXO 70 
employed is at frequency 4.096 MHz (or an harmonic thereof). The TCXO 
output may be used as the sample clock to the A/D converter 44 (or comparator); 
this acts to translate the signal to 1.028 MHz. This frequency is the difference 
between the 1 1th harmonic of 4.096 MHz and the input IF frequency 44.028 
MHz. The resulting 1.028 MHz IF is nearly one-fourth the sample rate, which is 
known to be nearly ideal in minimizing sampling type distortions. As compared to 
the I/Q sampling of Figure 2 A, this single sampler provides one channel of data 
rather than two, but at twice the rate. In addition, the data is effectively at an IF of 
1 .028 MHz. I/Q frequency conversion to near 0 MHz would then be implemented 
by digital means in the following processing to be described. The apparatus of 
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Figure 2A and 2B are competitive in cost and complexity; often component 
availability dictates the preferred approach. It will be apparent to those skilled in 
the an, however, that other receiver configurations could be used to achieve 
similar results. 

In order to simplify the following discussion, the following assumes that 
the I/Q sampling of Figure 2A is employed and that the snapshot memory 46 
contains two channels of digitized data at 2.048 MHz. 

Details of the signal processing performed in the DSP 32 may be 
understood with the aid of the flow chart of Figure 3 and the pictorial of Figures 
4A, 4B, 4C, 4D and 4E. It will be apparent to those skilled in the an that the 
machine code, or other suitable code, for performing the signal processing to be 
described is stored in EPROM 34. Other non-volatile storage devices could also 
be used. The objective of the processing is to determine the timing of the received 
waveform with respect to a locally generated waveform. Furthermore, in order to 
achieve high sensitivity, a very long portion of such a waveform, typically 100 
milliseconds to 1 second, is processed. 

In order to understand the processing, one first notes that each received 
GPS signal (C/A mode) is constructed from a high rate (1 MHz) repetitive 
pseudorandom (PN) pattern of 1023 symbols, commonly called "chips." These 
"chips" resemble the waveform shown in Figure 4A. Further imposed on this 
pattern is low rate data, transmitted from the satellite at 50 baud. All of this data is 
received at a very low signal-to-noise ratio as measured in a 2 MHz bandwidth. If 
the earner frequency and all data rates were known to great precision, and no data 
were present, then the signal-to-noise ratio could be gready improved, and the data 
greatly reduced, by adding to one another successive frames. For example, there 
are 1000 PN frames over a period of 1 second. The first such frame could be 
coherently added to the next frame, the result added to the third frame, etc. The 
result would be a signal having a duration of 1023 chips. The phasing of this 
sequence could then be compared to a local reference sequence to determine the 
relative timing between the two, thus establishing the so-called pseudorange. 

The above process must be carried out separately for each satellite in view 
from the same set of stored received data in the snapshot memory 46, since, in 
general, the GPS signals from different satellites have different Doppler 
frequencies and the PN patterns differ from one another. 
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The above process is made difficult by the fact that the carrier frequency 
may be unknown by in excess of 5 kHz due to signal Doppler uncertainty and by 
an additional amount due to receiver local oscillator uncertainty. These Doppler 
uncertainties are removed in one embodiment of the present invention by 
transmission of such information from a basestation 10 which simultaneously 
monitors all GPS signals from in view satellites. Thus, Doppler search is avoided 
at the remote 20. The local oscillator uncertainty is also greatly reduced (to 
perhaps 50 Hz) by the AFC operation performed using the base to remote 
communication signal, as illustrated in Figure 6. 

The presence of 50 baud data superimposed on the GPS signal still limits 
the coherent summation of PN frames beyond a period of 20 msec. That is, at 
most 20 frames may be coherendy added before data sign inversions prevent 
further processing gain. Additional processing gain may be achieved through 
matched filtering and summation of the magnitudes (or squares of magnitudes) of 
the frames, as detailed in the following paragraphs. 

The flow chart of Figure 3 begins at step 100 with a command from the 
basestation 10 to initialize a GPS processing operation (termed a "Fix Command" 
in Fig. 3) . This command includes sending, over a communication link 16, the 
Doppler shifts for each satellite in view and an identification of those satellites. At 
step 102, the remote unit 20 computes its local oscillator drift by frequency locking 
to the signal transmitted from the basestation 10. An alternative would be to utilize 
a very good quality temperature compensated crystal oscillator in the remote unit. 
For example, digitally controlled TCXOs, so-called DCXOs, currently can achieve 
accuracy of about 0. 1 parts per million, or an error of about 150 Hz for the LI 
GPS signal. 

At step 104 the remote unit's microprocessor 26 turns on power to the 
receiver front end 42, Analog to Digital Conveners 44 and digital snapshot 
memory 46, and collects a snapshot of data of duration K PN frames of the C/A 
code, where K is typically 100 to 1000 (corresponding to 100 msec to 1 second 
time duration). When a sufficient amount of data has been collected, 
microprocessor 26 turns off the RF to IF converter 42 and the A/D converters 44. 

The pseudorange of each satellite is computed in turn as follows. First, at 
step 106 for the given GPS satellite signal to be processed, the corresponding 
pseudorandom code (PN) is retrieved from EPROM 34. As discussed shortly, the 
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preferred PN storage format is actually the Fourier transform of this PN code, 
sampled at a rate of 2048 samples per the 1023 PN bits. 

The data in snapshot memory 46 is processed in blocks of N consecutive 
PN frames, that is blocks of 2048N complex samples (N is an integer typically in 
the range 5 to 10). Similar operations are performed on each block as shown in 
the bottom loop (steps 108-124) of Figure 3. That is, this loop is performed a 
total of K/N times for each GPS signal to be processed. 

At step 108 the 2048N data words of the block are multiplied by a complex 
exponential that removes the effects of Doppler on the signal carrier, as well as the 
effects of drifting of the receiver local oscillator. To illustrate, suppose the 
Doppler frequency transmitted from the basestation 10 plus local oscillator offsets 
corresponded to f e Hz. Then the premultiplication of the data would take the form 

of the function e-j2af c nT n= [ 0 , 1,2 2048N -1 ] + (B-l) x 2048N, where 

T= 1/2.048 MHz is the sampling period, and the block number B ranges from 1 to 
K/N. 

Next, at step 1 10, the adjacent groups of N (typically 10) frames of data 
within the block are coherently added to one another. That is, samples 0, 2048, 
4096, ... 2048(N-1) -1 are added together, then 1, 2049, 4097, ... 2048(N-1) are 
added together, etc. At this point the block contains only 2048 complex samples. 
An example of the waveform produced by such a summing operation is illustrated 
in Figure 4B for the case of 4 PN frames. This summing operation may be 
considered a preprocessing operation which precedes the fast convolution 
operations. 

Next, at steps 1 12-1 18, each of the averaged frames undergoes a matched 
filtering operation, whose purpose is to determine the relative timing between the 
received PN code contained within the block of data and a locally generated PN 
reference signal. Simultaneously, the effects of Doppler on the sampling times is 
also compensated for. These operations are greatly speeded, in one embodiment, 
by the use of fast convolution operations such as Fast Fourier Transform 
algorithms used in a manner to perform circular convolution, as presently 
described. 

In order to simplify discussion, the above mentioned Doppler 
compensation is initially neglected. 
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The basic operation to be performed is a comparison of the data in the 
block being processed (2048 complex samples) to a similar reference PN block 
stored locally. The comparison is actually done by (complex) multiplying each 
element of the data block by the corresponding element of the reference and 
summing the results. This comparison is termed a "correlation." However, an 
individual correlation is only done for one particular starting time of the data block, 
whereas there are 2048 possible positions that might provide a better match. The 
set of all correlation operations for all possible starting positions is termed a 
"matched filtering" operation. The full matched filtering operation is required in a 
preferred embodiment. 

The other times of the PN block can be tested by circularly shifting the PN 
reference and reperforming the same operation. That is, if the PN code is denoted 
p(0) p(l) ... p(2047), then a circular shift by one sample is p(l) p(2) .... p(2047) 
p(0). This modified sequence tests to determine if the data block contains a PN 
signal beginning with sample p(l). Similarly the data block may begin with 
samples p(2), p(3), etc., and each may be tested by circularly shifting the reference 
PN and reperforming the tests. It should be apparent that a complete set of tests 
would require 2048 x 2048 = 4,194,304 operations, each requiring a complex 
multiplication and addition. 

A more efficient, mathematically equivalent method may be employed, 
utilizing the Fast Fourier Transform (FFT), which only requires approximately 12 
x 2048 complex multiplications and twice the number of additions. In this 
method, the FFT is taken for the data block, at step 1 12, and for the PN block. 
The FFT of the data block is multiplied by the complex conjugate of the FFT of the 
reference, at step 1 14, and the results are inverse Fourier transformed at step 1 1 8. 
The resulting data so gotten is of length 2048 and contains the set of correlations 
of the data block and the PN block for all possible positions. Each forward or 
inverse FFT operation requires P/2 log2 P operations, where P is the size of the 
data being transformed (assuming a radix-2 FFT algorithm is employed). For the 
case of interest, B=2048, so that each FFT requires 1 1 x 1024 complex 
multiplications. However, if the FFT of the PN sequence is prestored in EPROM 
34, as in a preferred embodiment, then its FFT need not be computed during the 
filtering process. The total number of complex multiplies for the forward FFT, 
inverse FFT and the product of the FFTs is thus (2 x 1 1 + 2) x 1024 = 24576, 
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which is a savings of a factor of 171 over direct correlation. Figure 4C illustrates 
the waveform produced by this matched filtering operation. 

The preferred method of the current invention utilizes a sample rate such 
that 2048 samples of data were taken over the PN period of 1023 chips. This 
allows the use of FFT algorithms of length 2048. It is known that FFT algorithms 
that are a power of 2, or 4, are normally much more efficient than those of other 
sizes (and 2048 = 2* !). Hence the sampling rate so chosen significantly improves 
the processing speed. It is preferable that the number of samples of the FFT equal 
the number of samples for one PN frame so that proper circular convolution may 
be achieved. That is, this condition allows the test of the data block against all 
circularly shifted versions of the PN code, as discussed above. A set of alternative 
methods, known in the an as "overlap save 1 ' or "overlap add" convolution may be 
utilized if the FFT size is chosen to span a number of samples different from that 
of one PN frame length. These approaches require approximately twice the 
number of computations as described above for the preferred implementation. 

It should be apparent to one skilled in the art how the above process may 
be modified by utilizing a variety of FFT algorithms of varying sizes together with 
a variety of sample rates to provide fast convolution operations. In addition, a set 
of fast convolution algorithms exist which also have the property that the number 
of computations required are proportional to B log2B rather than B2 as is required 
in straightforward correlation. Many of these algorithms are enumerated in 
standard references, for example, H.J. Nussbaumer, "Fast Fourier Transform and 
Convolution Algorithms," New York, Springer- Verlag, C1982. Important 
examples of such algorithms are the Agarwal-Cooley Algorithm, the split nesting 
algorithm, recursive polynomial nesting algorithm, and the Winograd-Fourier 
algorithm, the first three of which are used to perform convolution and the latter 
used to perform a Fourier transform. These algorithms may be employed in 
substitution of the preferred method presented above. 

The method of time Doppler compensation employed at step 1 16 is now 
explained. In the preferred implementation, the sample rate utilized may not 
correspond exactly to 2048 samples per PN frame due to Doppler effects on the 
received GPS signal as well as local oscillator instabilities.- For example, it is 
known that the Doppler shift can contribute a delay error of ±2700 nsec/sec. In 
order to compensate for this effect, the blocks of data processed in the above 
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description need to be time shifted to compensate for this error. As an example, if 
the block size processed corresponds to 5 PN frames (5 msec), then the time shift 
from one block to another could be as much as ±13.5 nsec. Smaller time shifts 
result from local oscillator instability. These shifts may be compensated for by 
time shifting the successive blocks of data by multiples of the time shift required 
by a single block. That is, if the Doppler time shift per block is d, then the blocks 
are time shifted by nd, n=0, 1, 2, .... 

In general these time shifts are fractions of a sample. Performing these 
operations directly using digital signal processing methods involves the use of 
nonintegral signal interpolation methods and results in a high computation burden. 
An alternative approach, that is a preferred method of the present invention, is to 
incorporate the processing within the fast Fourier transform functions. It is well- 
known that a time shift of d seconds is equivalent to multiplying the Fourier 
Transform of a function by e~j2*fd, where f is the frequency variable. Thus, the 
time shift may be accomplished by multiplying the FFT of the data block by 
e -j27tndyT f for n=s0 , 1, 2, 1023 and by e -j2*(n-2048)d/T f for n= i024, 1025, 
... , 2047, where Tf is the PN frame duration (1 millisecond). This compensation 
adds only about 8% to the processing time associated with the FFT processing. 
The compensation is broken into two halves in order to guarantee continuity of 
phase compensation across 0 Hz. 

After the matched filtering operation is complete, the magnitudes, or 
magnitudes-squared, of the complex numbers of the block are computed at step 
120. Either choice will work nearly as well. This operation removes effects of 50 
Hz data phase reversals (as shown in Figure 4D) and low frequency carrier errors 
that remain. The block of 2048 samples is then added to the sum of the previous 
blocks processed at step 122. Step 122 may be considered a post processing 
operation which follows the fast convolution operation provided by steps 112- 
1 18. This continues until all K/N blocks are processed, as shown by the decision 
block at step 124, at which time there remains one block of 2048 samples, from 
which a pseudorange is calculated. Figure 4E illustrates the resulting waveform 
after the summing operation. 

Pseudorange determination occurs at step 126. A peak is searched for 
above a locally computed noise level. If such a peak is found, its time of 
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occurrence relative to the beginning of the block represents the pseudorange 
associated with the particular PN code and the associated GPS satellite. 

An interpolation routine is utilized at step 126 to find the location of the 
peak to an accuracy much greater than that associated with the sample rate (2.048 
MHz). The interpolation routine depends upon the prior bandpass filtering used in 
the RF/IF portion of the remote receiver 20. A good quality filter will result in a 
peak having a nearly triangular shape with the width of the base equal to 4 
samples. Under this condition, following subtraction of an average amplitude (to 
remove a DC baseline), the largest two amplitudes may be used to determine the 
peak position more precisely. Suppose these amplitudes are denoted Ap and 
Ap+i , where A p >A p+ i, without loss of generality, and p is the index of the peak 
amplitude. Then the position of the peak relative to that corresponding to Ap may 
be provided by the formula: peak location = p+ A p /(A p +A p+ i ). For example if 
Ap = Ap+ 1, then the peak location is found to be p+ 0.5, that is, halfway between 
the indices of the two samples. In some situations the bandpass filtering may 
round the peak and a three point polynomial interpolation may be more suitable. 

In the preceding processing, a local noise reference used in thresholding, 
may be computed by averaging all the data in the final averaged block, after 
removing the several largest such peaks. 

Once the pseudorange is found, the processing continues at step 128 in a 
similar manner for the next satellite in view, unless all such satellites have been 
processed. Upon completion of the processing for all such satellites, the process 
continues at step 130 where the pseudorange data is transmitted to the basestation 
10 over a communication link 16, where the final position calculation of the remote 
is performed (assuming Method 3 is utilized). Finally, at step 132, the majority of 
the circuitry of the remote 20 is placed in a low power state, awaiting a new 
command to perform another positioning operation . 

A summary of the signal processing described above and shown in Figure 
3 will now be provided. The GPS signals from one or more in view GPS 
satellites are received at the remote GPS unit using an antenna on the remote GPS 
unit. These signals are digitized and stored in a buffer in the remote GPS unit. 
After storing these signals, a processor performs preprocessing, fast convolution 
processing, and post processing operations. These processing operations involve: 
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a) breaking the stored data into a series of contiguous blocks whose 
durations are equal to a multiple of the frame period of the pseudorandom (PN) 
codes contained within the GPS signals. 

b) for each block performing a preprocessing step which creates a 
compressed block of data with length equal to the duration of a pseudorandom 
code period by coherently adding together successive subblocks of data, the 
subblocks having a duration equal to one PN frame; this addition step will mean 
that the corresponding sample numbers of each of the subblocks are added to one 
another. 

c) for each compressed block, performing a matched filtering 
operation, which utilizes fast convolution techniques, to determine the relative 
timing between the received PN code contained within the block of data and a 
locally generated PN reference signal (e.g. the pseudorandom sequence of the 
GPS satellite being processed). 

d) determining a pseudorange by performing a magnitude-squared 
operation on the products created from said matched filtering operation and post 
processing this by combining the magnitude-squared data for all blocks into a 
single block of data by adding together the blocks of magnitude- squared data to 
produce a peak. 

and e) finding the location of the peak of said single block of data to high 
precision using digital interpolation methods, where the location is the distance 
from the beginning of the data block to the said peak, and the location represents a 
pseudorange to a GPS satellite corresponding to the pseudorandom sequence 
being processed. 

Typically, the fast convolution technique used in processing the buffered 
GPS signals is a Fast Fourier Transform (FFT) and the result of the convolution is 
produced by computing the product of the forward transform of the compressed 
block and a prestored representation of the forward transform of the 
pseudorandom sequence to produce a first result and then performing an inverse 
transformation of the first result to recover the result. Also, the effects the Doppler 
induced time delays and local oscillator induced time errors are compensated for on 
each compressed block of data by inserting between the forward and inverse Fast 
Fourier Transform operations, the multiplication of the forward FFT of the 
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compressed blocks by a complex exponential whose phase versus sample number 
is adjusted to correspond to the delay compensation required for the block. 

In the foregoing embodiment the processing of GPS signals from each 
satellite occurs sequentially over time, rather than in parallel. In an alternative 
embodiment, the GPS signals from all in view satellites may be processed together 
in a parallel fashion in time. 

It is assumed here that the basestation 10 has a common view of all 
satellites of interest and that it is sufficiently close in range to remote unit 20 in 
order to avoid ambiguities associated with the repetition period of the C/A PN 
code. A range of 90 miles will satisfy this criteria. The basestation 10 is also 
assumed to have a GPS receiver and a good geographical location such that all 
satellites in view are continuously tracked to high precision. 

While several described embodiments of the basestation 10 show the use 
of a data processing component, such as a computer at the basestation in order to 
compute position information such as a latitude and a longitude for the mobile 
GPS unit, it will be appreciated that each basestation 10 may merely relay the 
information received, such as pseudoranges from a mobile GPS unit, to a central 
location or several central locations which actually perform the computation of 
latitude and longitude. In this manner the cost and complexity of these relaying 
basestations may be reduced by eliminating a data processing unit and its 
associated components from each relaying basestation. A central location, would 
include receivers (e.g. telecommunication receivers) and a data processing unit and 
associated components. Moreover, in certain embodiments, the basestation may 
be virtual in that it may be a satellite which transmits Doppler information to 
remote units, thereby emulating a basestation in a transmission cell. 

Figures 5 A and 5B show two embodiments of a basestation according to 
the present invention. In the basestation shown in Figure 5 A, a GPS receiver 501 
receives GPS signals through a GPS antenna 501a. The GPS receiver 501, which 
may be a conventional GPS receiver, provides a timed reference signal which 
typically is timed relative to GPS signals and also provides Doppler information 
relative to the satellites in view. This GPS receiver 501 is coupled to a disciplined 
local oscillator 505 which receives the time reference signal 510 and phase locks 
itself to this reference. This disciplined local oscillator 505 has an output which is 
provided to a modulator 506. The modulator 506 also receives Doppler data 



WO 97/14049 



PCT/US96/16161 



-26- 

information signals for each satellite in view of the GPS mobile unit and/or other 
satellite data information signals 511. The modulator 506 modulates the Doppler 
and/or other satellite data information onto the local oscillator signal received from 
the discipline local oscillator 505 in order to provide a modulated signal 513 to the 
transmitter 503. The transmitter 503 is coupled to the data processing unit 502 via 
interconnect 514 such that the data processing unit may control the operation of the 
transmitter 503 in order to cause the transmission of satellite data information, 
such as the Doppler information to a GPS mobile unit via the transmitter's antenna 
503a. In this manner, a GPS mobile unit may receive the Doppler information, the 
source of which is the GPS receiver 501 and may also receive a high precision 
local oscillator carrier signal which may be used to calibrate the local oscillator in 
the GPS mobile unit as shown in Figure 6. 

The basestation as shown in Figure 5A also includes a receiver 504 which 
is coupled to receive communication signals from the remote or GPS mobile unit 
via a communication antenna 504a. It will be appreciated that the antenna 504a 
may be the same antenna as the transmitter's antenna 503a such that a single 
antenna serves both the transmitter and the receiver in the conventional fashion. 
The receiver 504 is coupled to the data processing unit 502 which may be a 
conventional computer system. The processing unit 502 may also include an 
interconnect 5 12 to receive the Doppler and/or other satellite data information from 
the GPS receiver 51 1 . This information may be utilized in processing the 
pseudorange information or other information received from the mobile unit via the 
receiver 504. This data processing unit 502 is coupled to a display device 508, 
which may be a conventional CRT. The data processing unit 502 is also coupled 
to a mass storage device 507 which includes GIS (Geographical Information 
System) software (e.g. Atlas GIS from Strategic Mapping, Inc. of Santa Clara, 
California) which is used to display maps on the display 508. Using the display 
maps, the position of the mobile GPS unit may be indicated on the display relative 
to a displayed map. 

An alternative basestation shown in Figure 5B includes many of the same 
components shown in Figure 5A. However, rather than obtaining Doppler and/or 
other satellite data information from a GPS receiver, the basestation of Figure 5B 
includes a source of Doppler and/or other satellite data information 552 which is 
obtained from a telecommunication link or a radio link in a conventional matter. 
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This Doppler and/or satellite information is conveyed over an interconnect 553 to 
the modulator 506. The other input the modulator 506 shown in Figure 5B is the 
oscillator output signal from a reference quality local oscillator such as a cesium 
standard local oscillator. This reference local oscillator 551 provides a precision 
carrier frequency onto which is modulated the Doppler and/or other satellite data 
information which is then transmitted via transmitter 503 to the mobile GPS unit. 

Figure 6 shows an embodiment of a GPS mobile unit of the present 
invention which utilizes the precision carrier frequency signal received through the 
communication channel antenna 601 which is similar to the antenna 24 shown in 
Figure 1 A. The antenna 601 is coupled to the modem 602, which is similar to the 
modem 22 in Figure 1 A, and this modem 602 is coupled to an automatic frequency 
control circuit 603 which locks to the precision carrier frequency signal sent by the 
basestation described herein according to one embodiment of the present 
invention. The automatic frequency control circuit 603 provides an output 604, 
which is typically locked in frequency to the precision carrier frequency. This 
signal 604 is compared by the comparator 605 to the output of the GPS local 
oscillator 606, via interconnect 608. The result of the comparison performed by 
the comparator 605 is an error correction signal 610 which is provided to the 
frequency synthesizer 609. In this manner, the frequency synthesizer 609 
provides a higher quality, calibrated local oscillation signal over interconnect 612 
to the GPS down convener 614. It will be appreciated that the signal provided 
over interconnect 612 is similar to the local oscillator signal provided by 
interconnect 39 on Figure 1A to the converter 42; also, the converter 42 is similar 
to the GPS down convener 614 which is coupled to the GPS antenna 613 to 
receive GPS signals. In an alternative embodiment, the result of the comparison 
performed by comparator 605 may be output via interconnect 610a as an error 
correction to the DSP component 620 which is similar to the DSP chip 32 shown 
in Figure 1 A., In this instance, no error correction signal 610 will be provided to 
the frequency synthesizer 609. The automatic frequency control circuit may be 
implemented using a number of conventional techniques including a phase lock 
loop or a frequency lock loop or a block phase estimator. 

Figure 7 illustrates a particular sequence of power management according 
to one embodiment of the invention. It will be appreciated that there are numerous 
ways which are known in the art in order to reduce power. These include slowing 
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processed according to the method of the present invention, such as that shown in 
Figure 3. Then, after completing the processing of the GPS signal, the processing 
system is placed in a reduced power state as shown in step 23 (unless the 
processing system is also controlling power management as noted above). Then, 
in step 725 the communication transmitter in the modem 22 is returned to full 
power in order to transmit in step 727 the processed GPS signal back to the 
basestation 10. After completing transmission of the processed GPS signal, such 
as pseudorange information or latitude and longitude information, the 
communication transmitter is returned to reduced power state in 729 and the power 
management system waits for a delay of a period of time such as predetermined 
period of time in step 73 1. Following this delay the communication receiver in the 
modem 22 is returned to full power in order to determine whether a request is 
being sent from a basestation. 

Although the methods and apparatus of the present invention have been 
described with reference to GPS satellites, it will be appreciated that the teachings 
are equally applicable to positioning systems which utilize pseudolites or a 
combination of satellites and pseudolites. Pseudolites are ground based 
transmitters which broadcast a PN code (similar to a GPS signal) modulated on an 
L-band carrier signal, generally synchronized with GPS time. Each transmitter 
may be assigned a unique PN code so as to permit identification by a remote 
receiver. Pseudolites are useful in situations where GPS signals from an orbiting 
satellite might be unavailable, such as tunnels, mines, buildings or other enclosed 
areas. The term "satellite", as used herein, is intended to include pseudolite or 
equivalents of pseudolites, and the term GPS signals, as used herein, is intended 
to include GPS-like signals from pseudolites or equivalents of pseudolites. 

In the preceding discussion the invention has been described with reference 
to application upon the United States Global Positioning Satellite (GPS) system. 
It should evident, however, that these methods are equally applicable to similar 
satellite positioning systems, and in, particular, the Russian Glonass system. The 
Glonass system primarily differs from GPS system in that the emissions from 
different satellites are differentiated from one another by utilizing slightly different 
carrier frequencies, rather than utilizing different pseudorandom codes. In this 
situation substantially all the circuitry and algorithms described previously are 
applicable with the exception that when processing a new satellite's emission a 
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different exponential multiplier is used to preprocess the data. This operation may 
be combined with the Doppler correction operation of box 108 Figure 3, without 
requiring any additional processing operations. Only one PN code is required in 
this situation, thus eliminating block 106. The term "GPS" used herein includes 
such alternative satellite positioning systems, including the Russian Glonass 
system. 

Although Figures 1A, IB and 1C illustrate a multiplicity of logic blocks 
that process digital signals (e.g. 46, 32, 34, 26, 30, 28 in Figure 1 A), it should be 
appreciated that several or all of the these blocks may be integrated together onto a 
single integrated circuit, while still maintaining the programmable nature of the 
DSP portion of such a circuit. Such an implementation may be important for very 
low power and cost sensitive applications. 

It should also be appreciated that one or several of the operations of Figure 
3 may be performed by hardwired logic in order to increase the overall processing 
speed, while retaining the programmable nature of the DSP processor. For 
example, the Doppler correction capability of block 108 may be performed by 
dedicated hardware that may be placed between the digital snapshot memory 46 
and the DSP IC 32. All other software functions of Figure 3 may in such cases be 
performed by the DSP processor. Also, several DSPs may be used together in 
one remote unit to provide greater processing power. It will also be appreciated 
that it is possible to collect (sample) multiple sets of frames of GPS data signals 
and process each set as shown in Figure 3 while accounting for the time between 
the collection of each set of frames. 

A demonstration system, which is an example of an embodiment of the 
present invention, has been constructed that has verified the operation of the 
methods and algorithms described herein as well as showing the improved 
sensitivity possible by using these methods and algorithms. The demonstration 
system consisted of a GPS antenna and RF downconverter from GEC Plessey 
Semiconductors followed by a digitizer buffer board from Gage Applied Sciences, 
Inc. The antenna and downconverter perform the functions of 38, 40, 42, and 44 
of Figure 1 A and the digitizer buffer performs the functions 44, 46 and 48 of 
Figure 1 A. The signal processing was performed on an IBM PC compatible 
computer using a Pentium microprocessor, running under Windows 95 operating 
system. This emulated the functions of the DSP chip 32 and the memory 
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peripherals 34. Doppler information for satellites in view were provided to the 
signal processing software as inputs to the signal processing routines to emulate 
the functions of the modem and microprocessor 22, 24, 25, 26. 

The algorithms for this demonstration system were developed using the 
MATLAB programming language. A large number of tests were performed on 
live GPS signals obtained in various situations of blockage. These tests have 
verified that the sensitivity performance of the demonstration system was 
substantially superior to that of several commercial GPS receivers, that were tested 
at the same time. Appendix A provides a detailed listing of the MATLAB machine 
code that was used in these tests and is an example of the fast convolution 
operations of the present invention (e.g. Figure 3). 

In the foregoing specification, the invention has been described 
with reference to specific exemplary embodiments thereof. It will, 
however, be evident that various modifications and changes may be made 
thereto without departing from the broader spirit and scope of the invention 
as set forth in the appended claims. The specification and drawings are, 
accordingly, to be regarded in an illustrative rather than a restrictive sense. 
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APPENDIX A 



function [ sv , prange , snria , 6 n rout , evdoppler , outdAta ) «= 

gps< filename, s rate, code s Soppier, ao_pred, sample range) 

^function 

% [ sv , prange , snrin, snrout , svdoppler , outdata ] =gps ( filename , srate , codes , doppler , 

%no_pred, sampler ange ) 

% 

% N. Krasner Feb. 20 , 1996. - 

% CPrecision Tracking, Inc., San Jose, CA. 95117 
% 

%This function processes digitized data and provides pseudoranges as outputs 

%together with other statistical quanitites, as described below. 

t 

%Data to be processed is stored in a file named 'filename*. 
% 

%This function worXs on data from Gage digitizer at either sample rate 

% 3X2. 04 8 MHz (if erate-*0 ) or at sample rate 6X2.048 MHz (if srate-«l). 

%where it is assumed that IF is 35.42 -< 1400/45 )-4. 30B889 Hfiz. 

%The following are definitions of inputs and outputs of the routine: 

X 

ICodes are the SV numbers; if codes^O, then all codes will be used; 
% 

*Doppler « (dopper(l) doppler(2) ... doppler(n) span], specifies the mean 
%doppler of each of the n satellite plus a span to be seaxched over; if span 
%is missing or*0), then search is only done at the specified dopplers. 
%An algorithm is used to select the quantization size of doppler steps. 
%In some cases even if span is nonzero, only one doppler frequency per 
satellite %will be used. The search range is xspan relative to each doppler. 
% 

%no_pred is number of predetection frames 
% 

lAll data in the file is processed if samplerange is deleted. However, if 
%it is present it is two value and has the notation: 
% sampler ange={ lower_index, upper_index] ; 

%If these indices are not within range of that provided in the gage file, 

%then an error message is provided. 

% 

t Defaults are: if doppler has only one value, then doppler_span will be zero 
% if no_postd is not specified, then all data will be used 

% if no_pred is not specified, then pred will be set to 9 

% if no_pred is specified, then no_postd must be specified 

% 

I Outputs are codenumbers, pseudrangea , stirs and best doppler 

Iwith the pseudoranges expressed in nsec and snrs in dB. 

%A signal peak, must exceed IS dB threshold for detection. An 

% interpolation algorithm iB used to determine signal peak location. 

II f no codes axe found, the the returned function will be all zeros. 

% 

lOutdata is a vector of size 2048 that represents the output signal plus 
Inoise after postdetection integration for each detected satellite vehicle, 
%for the best doppler. This data was used to find pseudoranges 



%%%%% 

%First open the file and find the buffersize and indices to 



%%%% 
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tproceu, which may be all data if not specified i a input arguments 
f filevals, rate ]*gageopen( filename) ; 

buffs i2 e=fiieval3(3i; Itotal number of samples in file 

if nargm»=6, 

lower_index-samplerange< 1 ) ; upper_index-sampleranqe ( 2 ) - 
if upper_index>buffsize-l, 

error ( 'Sampler ange is out of bounds')* 

end 
else, 

lover_index=0;upper_index=buffsi2e-l ; 



end 
%%%% 



JeaTlf^H nUmbCr ° £ K d °PP ler bins ; the center of these bins are 

%called dopps and they are referenced to the center of the input dopplers 

if length (doppler )=* length < codes); span-=0; 
elseif length (doppler)— length ( codes ) + l # 

span«doppler { length (codes ) +1 ) ; 
else, 

end err ° rrLen ? th of <*<>PPler vector must equal length of codes, or 1 more • ) ; 

if span<0, error('Span must be greater than zero');end 

if srate«0, s_per_ frame- 6 ; else, s_per_f rame«12 ;end 

if nargin«=3, no_pred-min (9, fix (buff size/ ( 1024* s_per frame))) -end 
no_postd=fix(buffsiie/(no_pred*1024*s _per_f rame) ) ; ~ 

if no_pred<l, error ( 'Number of predetection frames must be at least iM-end 
no_j>ostd<l, error( -Number of postdetection frames must be at least i : );end. 

deltaf-(l/3)*1000/no_pred; %yields about 1.5 dS loss which is made up 

%worst case by two independent detections for 
%when true freq is midway between bins 

if span<delta£, 

dopp«=0; %no need to quantize doppler 
else, 

dopps=[0:-deltaf :-span] ; 

dopps«sort{ [dopps deltaf : del taf: span] ) ; 

end 

%«% 

_ — — — _ % % % 

%Now get the FTT of codes and put them in a matrix of size lengths codes )X2 04 8 

if codes~0, codes«l:28;end %«earch all codes in this case 

cod«array«zero8< length (codes), 2 046); 
for k« 1 : length < codes ) , 

codear*ray(k, : )■ 9P s_rea amp led ( codes < Jc) ) ; 

end 

III ~ . _ %%% 

%In the following we input data in blocks equal to a number of PN 
tframea equal to nojred. For example if no_pred-7, the the number 
lot samples from the gage is 2048«3-7«43008, bo that in this case 
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*a full 4 Mbyte collection would cake .bout 93 ite -,-^. 
*a«r« tl oa then performs the orec suraa-or f" ° Ch 

Squared operation and station ^trlTlou "™-"°^ el «iM. Mg - 
*an= all satellite vehicles specified/ CkS f ° r 811 d »??i«-'s 

biocksize=1024's_per_irameT.o _pred; 

■start=lower_index-blocksi2e ; 
lc-length( codes ) ; ld-=lencth ( dopps ) ; 
prriata-zerosf le*ld,2048) ; 

no_blockG-f ix( (buf f size-10-s nc- < r ,.„, /u , . 

««jnples s _Pe-_-r«ne)/Mocksize); » extra 6Q or 120 

» % 60 r or V ?50 aP "" Ve /, U ° WS i " pui8e "«PO«-e to 
for n-l:no_block S , * 120; 9e ° lf -2-*>«e m-fii e 

start=st«rt+blocksi ze ; 

if till STk ' " , ?* Sr ' ad ' f ileV * ls ' - blocksire + 10 • s _per frame , • 

"■ i ' J>se,dat4 ' s " te,; t rsrs" data to - 

r^^f 1:2 ° +204 n^:e d, o V e- % sv 5e " — — — — — 

for d=l:ld, %Iter«te over dopplers ==°e<k) 

df=doppler(k)*dopps(d) ; 

rdata=data . -exp ( - 3 -2 .pi«df • f 0 : 2048-no_pred-l ) 17 .048eS J • 
rdata-sun(reshap«(rdata.204B,no ored) 7 )- tdo .. Zl * 

%corapeneate fcr tune slip per block which * * f ■ 
%noting that fdoop/fo is^pler tL " ° ^ ^ 

runesl i p=le-3-no_j ) red*df/1575.42e6- second 

P||"e-exp(-3-2-pi-tiiiieslip-(n-l)-[ 0: 1024]-le3)- 
rdata(lsl025)-rdata( 1:1025). -phase- 

SUISSS^Ltt 1- " 1 10215 S2MB 1 - -~3n*-«.f io2* aa , i , 

rdata-real(rdata).-2*ijnag(rdata| .-2: *matmitud» 
prdata<, k -l ) .ld*d. : )-prd.%., ( k-i ) .i; +d .:, + rd.^r^ 6 

end 

end 

end 



%%% 



*Now we look for peaks exceeding threshold 
for k-l:lc. %iook at: all SVe and dopps 

^£iiEi.;r il,ld, ' s,j %data for sv#k - ■» ^» 

%yields Pfa less than le-4 per using all 

< - . A %29 av's and 10 dopplere dc- sv 

mdx*find(2<threshold); pp ^ sv 

d-i!r;? rt<meen<mean<zUndt> ^ 2 > ,,; defined noise getting rid of aicmal 
d-z.O^.n.aximrotatet z.-l ) ,mrotate«z, 1 » , ; idetect.on 2ft IlL^natxng^ 

f _ Jw j - , , ftpicicet fence problem 

[™xd.i„dd)=n,ax ( d-,; %£indg max and loe « ioB for ^ d<J ^ ler 
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[maxsv,indsv)*max(maxd); %find global peak for this sv 

if maxsv> threshold, %then we detect this sv! 

r=(maxsv/noiser2; %peak to noise away from peak 
smn=sqrt ( r/( no_pred* 2*no_postd) } ; %anrin via N.K analysis 
snxin*[snrin 10*logl0( enin) ] ; - y 6 

%now find anrout at peak 
S nout-(no - pred-2-no_postd ) *snin-2/(l + 2-no_pred« S nin); ^i^N K anal 
snxout={snrout 10«logl0< snout )] ; * nai * 

sv»[sv codes (k) ] ; 

d=z<indsv, : ); %data for beat doppler 

outdata=(outdata d); 

svdoppler- ( svdoppler doppler ( k > +dopps ( indsv ) J ; %beat doppler 
indO*indd(andsv) ; %best sample for best 

indsO«(indO-4:ind0+4); %block of 9 data samples about indO- 
nxndssf ind{ inds0<*0 ) ; 
pinds=£ind(indsO>204B) ; 

if length<ninds)>0, %keep data to interpolate within range 

ind«0( ninds)=inds0(ninds)+2048; 
elaeif length ( pinds >>0, 

indsO ( pinds }«inds0( pinds ) -2046 ; 

end 

dint«interp(d(inds0) ,10) ; %go to rate 20 Mflz 
dint=interp(dint,10); % go to rate 200 MHz 

dint-interp(dint,5) ; %go to rate 1 GBz 

(amax,imax) -max (dint); %find max 

in<taiax=min( find (dint> 0.8 *araax) ) ; %max above 80% 
tmax==indO+(indmax-2001)/500; %oseudorange in samples 
if tmax>*2048, 

tmax»tmax-204 8 ; 
eleeif traax<0, 

tma x« una x+ 2 0 4 8 ; 

end 

tma.x«tjoax*500; %p«eudorange in units of nsec 
pr ange= I pr ange tmax J ; 

end 

end 

if length <av)«*o, %then there were no detections 

sv-0 ; praage*0 ; snrin-0 ; snxout«0 ; evdoppler=0 ; outdata«0 * 
end ' 

f close ( 'all ) ; 



function Y-gps codes (num) 

ly^gpacodea ( num ) 

% 

*% H. Kraaner Jan. 22, 1996. 

% OPrecision Tracking, Inc., San Jooe, CA. 95117 

%Returns all 1023 bit codes associated with satellite in an 
%array y of size 37X1023. Each row is the code associated 
%with the satellite having the corresponding ID number. That 
%is row 1 is the pn (Gold) code of satellite 1, etc. 

%If num is specified then only one code with that number is 
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% returned 

yl=pngen( [3, 10] , 102 3, ones ( 1, 10) ) ; 
y2=pngen< [2, 3,6,8,9,10] , 1023,ones{ 1 , 1 0 > ) ; 

delays=[5,6 7,8,17,16,139, 140,141,251,252 254 255 r>5< ->c Q , 
delays=(delay S , 469,470,471,472,473,47 , 0^5 2 513 SU^sfi; 
^elays-fdelays, 859 860 861 862 663 950 947 948 950J; ^ 

if nargin-«0, 

for k»l:iength(delays), 

y(k, : > =xor(yl, rotate (y2, delays* X) ) ); 

end 
else, 

y-xor <yl , rotate < y2 , delays ( nmn) ) ) ; 
end 



function <-gps_rsnp( code number) 
* N. Xrasner Feb. 16, 1996. 

% ©Precision Tracking, Inc., San Jose, CA. 95117 

^e Ura ! ^ resam P led ^urier transform of specified 1023 bit 

%GPS code, resampled to at rate 2048 JcBz and weighted 

%with a -sane- function to emulate the transform of the 

^sequences consisting of a set of square pulses, which 

lis ideally bandpass filtered to the first nulls (±1023 Jtflz). 

%of "ze n ^048 8 °' then aU W ' COdeS *" rCtUrned in * Mtei * 

% 

^-^h t ^ n °u° r l!! ed d * ta contains real «d imaginary components 
»e«ch of which 1S rounded to 4 bit signed integers in the range 
M-7,7]. This is the form stored in EPROM . This quantisation 
SSLI eqUired B *" 0ra9e ' bUt r " Ult9 " "Eligible loss of 

% 

%As in the function gpscodes, the returned value contains 37 
*rows each corresponding to the corresponding GPS code (row 1 

^rien^h^Ifl 1 ' if * Xl COd ° 8 M "^"^ci.Each row is 

*!lJ™r * nd/ as stated contains complex words whose 

*camporeents are integers in the range [-15,15], 

%Note that this function calls the function gpscodes to create 
%the various Gold codes. 

if codenumber-«0, 

y^gpscodes ; 
else, 

y*» gpscodes ( codenumber ) ; 

end 

% returns matrix 37X102 3 of GPS codes 

y " 2 * y " 1; %convert.to *1 

if codenumber»o, 
z=2eros( 37,2048); 
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for k=l;37, 

iyl=fft< stretch (y<fc, : ) ,2 ) ) ; 

*yll = [fyl(l:1024) .*sinc< [0:1023)/1023) 0 0J- 
fyll-ffyll fyl( 1025:2046) .'sine ([ 1022:-! :1 J/1 023J J- 
2(k f r)-conj(fyll); %need conjugate for matched filter 

end 

else, 

fyl^fftf stretch(y,2) ); 

fyll = [fyl(l:1024) .•sinc([0;1023]/1023) 0 01- 
tyll-[fyll f yK 1025:2046). •sinc(fl022:-l:l ]/l023ll- 
2=conj(fyll) ; %need conjugate for matched *ilte- ' 
end 
end 



z=round(z/20) ; %since max is 171 
rz=min(real(2) ,7 ) ;rz«max(r2,-7) ; 
iz«min( imag( z ) , 7 ) ; iz«max< iz , -7 ) ; 
z^r2+j«iz. ; 



9 but few values are above 140 
%limit range 
%limit range 



function <Uta_out«if_2_base(data_in, rate) 

%data_out=if _2_base ( data_in) 

% N. Krasner Feb. 20, 1996. 

% CPrecision Tracking, Inc., San Jose, CA. 95117 
% 

%If rate is missing or equals 0, then the input rate is assumed 
%to be 3-2048 *fl2. If rate equals 1, then it is 6*2048 JcBz. 

%Converts data from GAGS digitizer sampled at either 3*2048 kBz 

lor 6*2048 kBz with IF 6.144 Mfl2 (3*2.048) centered at 

%35.42 -(1400/45)-4. 308889 MB2 to a complex <I,Q) stream 

*of data with sample rate 2.048 Mfiz centered at 0 

Ifrequency (nominally). This routine does not compensate for 

*any doppler shifts on SV's, nor any LO drifts. Note that the output 

%data length is either 1/3 or 2/3 that of the input. 

if nargin««l, rate«0;ead 

lfreq - 35.42 - (1400/45); %in MHz 

if rate— 0, %case where sample rate is 3*2048 xBz 

srate-3 -2.048; % in KHz % s ample rate for GAGE 
{ u , v ) -size ( data_in ) ; i f u> v ; dat a_in«da ta_in ■ ; end ; Id- length ( data in ) ; 
c«ejrp(-}*2*pi*(0:ld-l]*ifreq/erate); %LO ~ 
- data-out-data-in. *c; clear data in c %eonvert: to 0 

h-firl(50,.9/3.072); " tconvem to 0 

data_out-filter<h,l,data_out); % filter to =1 MHz 

data_out-data_out(l:3;ld); %decimate to yield 2.048 MB 2 srate 

* **' ^ * * „ %caae where sample rate is 6*2048 Jcflz 

srate«6*2.048; % in KHz 

[u,v)«size<data_in>;if u>v; data-in-data in* ; end ^d-lengthl data in)- 
c-exp(-3«2«pi*t0ild-l)*ifreq/srate); \Lo " 
data_ou*«data_in.*c; clear data in c %convert to o 

h-firl<100, .9/6.144); " convert to 0 

data_out-filter<h,l, data-out); tfilter to *l HBz 

data_out-data_out:(l:6:ld); Idecimate to yield 2.048 KBz srate 
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end 



function y«nrotate(z f B) 

%y«inrotate( x,m) ; matrix rotation by m units 
I 

X N. Krasner Feb. 20, 1996. 

% OPrccision Tracking, Inc., San Jose, CA. 95117 
m=round(m) ; 
Tn*rem(m, lx) ; 
if n>0, 

y=[x< : , lx-m+1 :lx) x( : , 1 : lx-m) ] ; 
elseif m<0, 

y*[x( : , l-m:lx) x(:,l:-m}]; 
else , 

y=x; 

end 



function y-pngen( tap* ,npts, fill) 

Ipngen ( taps , npts , f ill ) 

% 

% K. Rrasner Feb. 13, 1995. 

I ©Norman F. Krasner, San Carlos, CA. 94070 

% 

% PNGEN generates 'npts' number of points of the peuedo random 

\ sequence defined by a shift register of length n with feedback 

taps 

% "taps' and initial fill 'fill'. PN(npts, taps, 

fill) returns a 

% vector of length npts containing this sequence. 

* A typical example [1,6,8,14] for R14, filial 

ieros(l,13) ] 

% 

% If a single number is given in "taps" it is interpreted as 

% the octal value specifying the taps, as in Peterson and Weldon. 

For 

I example a code 2767 (octal) = 010 111 110 111. We delete the 
leading 0 

I to get the polynomial 10 111 110 111. The feedback taps for this 
case 

% are found by ignoring the first one and choosing the remaining 
" ones * 

% to get position. For this case the taps are thus 2 34 5689 10. 

if length ( taps )»« 1 , 

taps»«printf ( ' t . Of ' , taps ) ; ntaps* [ ] ; 

if taps( 1 )*»• 1 ' r ntaps- [1] ; 

elseif taps( 1)««*2 ' , ntaps- { 10]; 

elseif taps( 1 )«•= ' 3 ' , ntaps* (1 1 ]; 

elseif taps( 1 )«« ■ 4 * , ntaps* (1 0 0]; 

elseif taps( 1 )«* * 5 * , ntaps-fl 0 1]; 

elseif tapa( 1 )« ' 6 * , ntaps s (l 10); 

elseif taps( 1)«**7 ' , ntaps* [1 11]; 
end 
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for k=2 :length( taps) "39" 

if taps{ k)«= *0 ' , ntaps*=( ntaps 0 0 0] 

elseif taps(k)«='l ■ , ntaps= [ r.taos 0 0 1] 

elseif taps(k>«=*2 ' , ntaps=[ntaos 0 10] 

elseif taps(k>==3* , ntaps= [ ntaos Oil] 

elseif taps(k)»'4* r ntaps=(ntaps 10 0] 

elseif taps(k)«*5*, ntaps=[ntaos 10 1] 

elseif taps(k>«='6 ■ , ntaps=[ntaps 110] 

elseif taps(k)«= '7 1 , ntaDs*(ntaps 1 1 l i 

end 1 

end 

ntaps(l)=f ]; 
taps=f indfntaps ) ; 

end 



n«max< taps ) ; 

if nargin<=2, fill=[l zeros ( 1 ,n-l )] ;end 
if aargin-=l, npts« ( 2*n)-l ;end 
y - zeros (1, opts); 
x « zeros(l,n); 
x = fill (length (fill) :-l:l); 
for i=l:npts 
y(i) = x<n); 

x~[rem(3um(x(taps) ) ,2) , x(l:n-l)l; 

end 



function (file_Tals, aajople^rate ] - gageopen ( f ilena*e ) 

* Read gage formatted file, with particular number of samples 

* function [file_vals, sample_rate] « gageopen ( filename) 

* 3. Wilson Jan. 13, 1996. 

I OPrecision Tracking, Inc., San Joae, CA. 95117 
file__vala(l) = fopen< filename, ■r 1 ); 

(sample rate, sample_depth, head_vala, operation mode] - 
rdgghead(file_valB(l) ); 1 

file_valo(2j - operation mode; 
file_vals(3) «= sample_depth ; 
file_valo(4:6) • head_vals( 1 : 3 ) 

function fbuf, count] - gageread(f ile_Tal., start, nu.sa.ple8) 

* Read gage formatted file, with particular number of samples 

i B . WU^aan!^'; " '•^<"^-.-««.»-— Pi") 

I OPrecision Tracking, Inc., San Jose, CA. 95117 



% fil«s_vals 

* 1~- file id 

L^i V f**"^ 0 "-* 0 ** - if -1; file is memory image interleave 
sample_depth long; =2 memory image 
% 3 - aample_depth 
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\ 

ns2 = numsamples/2 ; 
fid = file_vils( 1 ) ; 
if (fJ.le_vals(2 ) == 2) 

buf = zeros { numsa.nples, 1 ) ; 

fresult = fseek<fid,512*start,-M - * m o V „ . 
data portion of file * 1 ' * move to "quested point in 

else (bUf ' C ° UntJ * fr€ad{fid ' numsaiaples, 'uchar'); 
buf » zeroa(ns2,2 ) ; 

fresult = fseek<fid,512+<3tart/2i -m< » 
point in data portion of file rt/2) ' * mov * *o requested 

[buf(:,l), clj = f read (fid, ns2, 'ucharM- 

fresult = faeekf fid, 512*< start/2 )+file vaia<3i i \ . 
BANK B, portion of memory ijnage ^ie_yais< 3 > , -1 , ; % move to 

[buf(:.2>, c2] = freadtfid, ns2, 'ucharM- 

count = cl + c2; 

buf = reshape<buf ■ , numsamples, 1 1 • 

end 
% 

buf - buf - 12 8; 



* Read gage formatted header 
% B. Wilson Jan. 13, 1996. 

% CPrecieion Tracking, inc., San Jose, CA. 95117 

% 
% 

t sample_rate = 

tl.2.4,5.1o7l2.S # 20,2S.30,40.50,60.100.I20.125.150,200.250j; 
fresult - fs«sk(tid,C,-l); % rewind the file 

^"P " fread(fid, 14, char - ); 

• » •et■tr(tlIlp• ) ; 

fresult - f«*»k< fid. 16.-1); % move to name variable 
tap - freed ( f id , 9 , • char • ) ; 
•nan* ■ aetotr (tjnp- ) ; 

varia^e" £ '~ k < fid < 287 <"1 > = * »°ve to index to a^ple rate table 

sanple_rate_index = f read( fid, 1 , • intl6 • ) - 
if <»aaple_rate_index «= 42) » 42 is external" 

»a«ple__rate « -1; 

else 



end 
% 



■ample_rate = aample_rate_table< aample_rate_index - 17); 



fresult . faeek(fid,2B9,-l); % move to index to operation mode variable 
operationmode - tread< fid. 1 , " intl6 • ) P«a*ion nuxie variable 

*% if (operation_mode « 1) 
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sa.xple__rate = sample_rate / 2; 

\\ end 



fresult = fseekffid, 301,-1); % move to sample depth variable 
sajnple_depth = f readf f id , 1 , • int32 ' ) - ae P- n va.iaile 

- - ( opera tion_mode « 1 ) 

sample_depth - sampie_depth / 2 % eeens like a "bug- to me 
if (sample_depth 4194304) 



sarople_depth = sample_depth /2 
end 

end; 



fresult = fseek(fid.313,-l>; % move to starting address 

head_vals = f readf f id, 3 , * int32 ' ) ■ 

% 

fresult - fseek( fid, 339,-1); % move to sample depth variable 
resolution_12_bits = f read< f id, 1 , • mtl 6 * ) ; ~ variable 

^f!? 1 * L fs€fiJc < £id '345,-l, ; % move to sample depth variable 
sampie_offset » f read ( f id , 1 , * intl 6 • ) ; " 

fl^V feee *' fid ' 349 - 1 >; * »°ve to sample depth variable 
sanple_bit6 * fread(fid,l, *intl6'); ~ 



function y*rotate (x ,q) 

%y*rot*te(x,q) 

%Rotates vector by q places 

% N. Krasner Nov. 28, 1995. 

% ©Norman Krasner, San Carlos, CA. 94070 

last q half deleted r0tates a Vector ^ ^""hanging the first and 

JpLced It C thf 'J** VC r° r " ° dd in len ^ h the ^ d Point is 

*piaced at the beginning of the new vector. 

%If q is present and >0 this rotates the vector x right by q positions; 

then^ 1 1 ° 1 01 ' thCn rotate < x ' 2 > Yields y-[l 0 110). If q is <0 
%this rotates x to the left. 

[n»,n)«size(x) ;mn«max(m,n) ; 
if nargin>l f q«r em (q,mn) ; end 

UrTiLt « rror ' ,Tai » Unction works only for vectors, not 

if m> 1 , x»x ' ; end 

if nargin«*l, 

if rem (mn, 2 ) ==0 , 

pivot- l+mn/2; 
else pivot *( ran+1 ) /2 ; 
end 

y* {*< pivot :mn) x( 1 :pivot-l ) J ; 
elseif aba<q)>0, 

if q<0, q-q^nn;end %fix for Bni£t8 left 

y«(x(mn-q+l:mn) x(l:mn-q)); 
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else y=x; 
end 

end 

if n> I , y*y ' ; end 

function z«sinc(x) 
%2*sinc<x) 

4 N . Krasner April 21, 1991 

% ONorman Krasner, San Carlos, CA. 94070 

I 

(Definition: sinc(x) = sin( pi*x) /(pi'x) , where x is a vector /matrix. 
%sinc(x)=l if x=0 

indl=f ind( x«0 ) ; 
ind2=f ind( x-=0 ) ; 
z«x; 

x( indl )«ones( 1 , length( indl ) ) ; 
x(md2)=sin(pi*x(ind2 ) ) . / (pi*x( ind2 ) ); 
z{:>=x; 



function z -st re t ch <■ ignal , period) 
(function z=stretch< signal, period) 
X 

% N. Krasner March 12, 1995 

t ONorman Krasner, San Carlos, CA- 94070 

(This function stretchs a signal oy replacing each element 
(with m identical elements; m is specified by "period.* 

[u,v]=size( signal) ; 

if u>l 4 v>l error ( 'This function works only for vector inputs.'); end 

if u>l, signal«signal ' ;end 

a«ones ( period , 1 ) * signal ; 

z-&( : ) * ;z-=conj ( z ) ; 

if u>l, z«*conj ( z ' ) ;end 



function y*xor(a,b); 

% function y»xor(a,b); 

(Exclusive or, term by term, of two vectors 

( N. Krasner June 13, 1994 

% ONorman Krasner, San Carlos, CA. 94070 



y~(afc-b ) | (bfc-a); 
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Claims 

1 . A GPS receiver apparatus, comprising: 

an antenna for receiving GPS signals at an RF frequency from in 

view satellites; 

a downconverter coupled to said antenna, said downconverter 
reducing the RF frequency of said received GPS signals to an intermediate 
frequency (IF); 

a digitizer coupled to said downconverter and receiving said IF 
GPS signals, said digitizer sampling said IF GPS signals at a predetermined rate to 
produce sampled IF GPS signals; 

a memory coupled to said digitizer, said memory storing the 
sampled IF GPS signals; and 

a digital signal processor (DSP) coupled to said memory, said DSP 
performing fast convolution. 

2 . A GPS receiver as in claim 1 further comprising a communication 
antenna and a receiver coupled to said communication antenna and to said DSP, 
said receiver for receiving a data signal containing satellite data information. 

3 . A GPS receiver as in claim 2 wherein said satellite data information 
comprises a Doppler information of a satellite in view of said GPS receiver. 

4. A GPS receiver as in claim 3 wherein said satellite data information 
comprises an identification of a plurality of satellites in view of said GPS receiver 
and a corresponding plurality of Doppler information for each satellite of said 
plurality of satellites in view of said GPS receiver. 

5 . A GPS receiver as in claim 2 wherein said satellite data information 
comprises data representative of ephemeris for a satellite. 

6. A GPS receiver as in claim 1 further comprising: 

a local oscillator coupled to said downconverter, said local 
oscillator providing a first reference signal. 

7 . A GPS receiver as in claim 2 further comprising a local oscillator 
coupled to said downconverter, said local oscillator providing a first reference 
signal and wherein said receiver receives a precision carrier frequency signal 
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which is used to calibrate said first reference signal from said local oscillator, 
wherein said local oscillator is used to acquire GPS signals. 

8 . A GPS receiver as in claim 3 wherein said DSP compensates said 
sampled IF GPS signals using said Doppler information and wherein said fast 
convolution operations provide a pseudorange information. 

9 . A GPS receiver as in claim 1 further comprising a power 
management circuit coupled to said downcovener and to said digitizer, wherein 
after said IF GPS signals are stored in said memory, said power management 
circuit reduces power consumed by said downconverter and said digitizer. 

10. A GPS receiver as in claim 8 further comprising a transmitter 
coupled to said DSP, said transmitter for transmitting said pseudorange 
information. 

11. A GPS receiver as in claim 2 further comprising a transmitter 
coupled to said DSP said transmitter for transmitting a latitude and longitude 
information. 

1 2. A method for using a GPS receiver, said method comprising: 
receiving GPS signals from in view satellites; 

digitizing said GPS signals at a predetermined rate to produce 
sampled GPS signals; 

storing said sampled GPS signals in a memory; 

processing said sampled GPS signals by performing fast 
convolution operations on said sampled GPS signals in said GPS receiver. 

13. A method as in claim 12 further comprising: 

receiving a data signal containing satellite data information. 

14. A method as in claim 13 wherein said satellite data information 
comprises a Doppler information of a satellite in view of said GPS receiver. 

15. A method as in claim 14 wherein said Doppler information is used 
to compensate said sampled GPS signal and wherein said processing further 
comprises preprocessing and postprocessing operations. 

16. A method as in claim 15 wherein said fast convolution operations 
provide a pseudorange information and wherein said processing further comprises 
preprocessing and post processing operations. 

17. A method as in claim 13 wherein said satellite data information 
comprises data representative of ephemeris for a satellite. 
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18. A method as in claim 17 wherein said fast convolution operations 
provide a pseudorange information and wherein said ephemeris and said 
pseudorange information are used to calculate a latitude and longitude of said GPS 
receiver, 

19. A method as in claim 18 wherein said latitude and longitude is 
displayed to a user of said GPS receiver. 

20. A method as in claim 18 wherein said latitude and longitude are 
transmitted by said GPS receiver. 

21 . A method as in claim 12 wherein said GPS signals originate from 
pseudolites. 

22. A method as in claim 12 wherein said GPS signals originate from 
orbiting satellites. 

23. A GPS receiver as in claim 1 wherein said GPS signals originate 
from pseudolites. 

24. A GPS receiver as in claim 1 wherein said GPS signals originate 
from orbiting satellites. 

25. A method of determining pseudoranges in a global positioning 
satellite (GPS) receiver, comprising: 

receiving GPS signals from one or more in view GPS satellites 
utilizing an antenna coupled to a downconverter, said GPS signals comprising 
pseudorandom sequences; 

buffering the received GPS signals in a digital snapshot memory; 

processing the buffered GPS signals for one or more of the in view 
GPS satellites in a digital signal processor by: 

breaking the buffered data into a series of contiguous blocks whose 
durations are equal to a multiple of the frame period of the pseudorandom (PN) 
codes contained within the GPS signals, 

for each block, creating a compressed block of data with length 
equal to the duration of a pseudorandom code period by adding together 
successive subblocks of data, said subblocks having duration equal to one PN 
frame, such that the corresponding sample numbers of each of the subblocks are 
added to one another; 

for each compressed block, performing a convolution of the 
compressed block's data against the pseudorandom sequence (PRS) of the GPS 
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satellite being processed, said convolution being performed using fast convolution 
algorithms, said convolution producing a result; 

performing a magnitude-squared operation on the results created 
from each of said convolutions to produce magnitude-squared data; 

combining said magnitude-squared data for all blocks into a single 
block of data by adding together such blocks of magnitude-squared data, such that 
the corresponding sample numbers of each of the magnitude-squares from the 
convolution are added to one another; and 

finding the location of the peak of said single block of data to high 
precision using digital interpolation methods, where the location is the distance 
from the beginning of the data block to the said peak, and the location represents a 
pseudorange to a GPS satellite corresponding to the PRS being processed. 

26. The method of claim 25 wherein the fast convolution algorithm 
used in processing the buffered GPS signals is a Fast Fourier Transform (FFT) 
and the result of the convolution is produced by computing the product of the 
forward transform of said compressed block and a prestored representation of the 
forward transform of the PRS to produce a first result and then performing an 
inverse transformation of said first result to recover said result. 

27 . The method of claim 25 wherein the fast convolution algorithm 
used in processing the buffered GPS signals is a Winograd algorithm. 

2 8 . The method of claim 26 wherein the effects of Doppler induced 
time delays and local oscillator induced time errors are compensated for on each 
compressed block of data by inserting between the forward and inverse Fast 
Fourier Transform operations, the multiplication of the forward FFT of said 
compressed blocks by a complex exponential whose phase versus sample number 
is adjusted to correspond to the delay compensation required for said block. 

29 . The method of claim 25 wherein the digital signal processor is a 
general purpose programmable digital signal processing chip which is executing 
stored instructions. 

30. The method of claim 25 wherein the fast convolution algorithm 
used in processing the buffered GPS signals is a Agarwal-Cooley algorithm. 

3 1 . The method of claim 25 wherein the fast convolution algorithm 
used in processing the buffered GPS signals is a split nesting algorithm. 
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32. The method of claim 25 wherein the fast convolution algorithm 
used in processing the buffered GPS signals is a recursive polynomial nesting 
algorithm. 

33. The method of claim 25 further including the step of: 
determining that said peak is valid by determining whether said peak 

exceeds a predetermined threshold. 

34 A tracking process utilizing global positioning system (GPS) 
satellites for determining the position of a remote sensor, the process comprising 
the steps of: 

receiving and storing GPS signals at said remote sensor from a 
plurality of in view GPS satellites; 

computing in the sensor, pseudoranges utilizing said GPS signals, " 
said computing comprising digital signal processing using fast convolution 
techniques on stored GPS signals; 

transmitting said pseudoranges from said sensor to a base station, 
said base station being provided with GPS satellite ephemeris data; and 

receiving said pseudoranges at said base station and utilizing said 
pseudoranges and said satellite ephemeris data to compute a geographic location 
for said sensor. 

35. The tracking process of claim 34 wherein the step of computing the 
pseudoranges further comprises: 

storing the received GPS signals in a memory; 

processing the stored GPS signals for one or more of the in view 
GPS satellites in a digital signal processor by: 

breaking the stored data into a series of contiguous blocks whose 
durations are equal to a multiple of the frame period of the pseudorandom (PN) 
codes contained within the GPS signals; 

for each block, creating a compressed block of data with length 
equal to the duration of a pseudorandom code period by coherently adding together 
successive subblocks of data, said subblocks having duration equal to one PN 
frame; 

for each compressed block, performing a matched filtering 
operation to determine the relative timing between the received PN code contained 
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within the block of data and a locally generated PN reference signal, said matched 
filtering operation utilizing said fast convolution techniques; 

determining said pseudorange by performing a magnitude-squared 
operation on the products created from said matched filtering operation and 
combining said magnitude-squared data for all blocks into a single block of data by 
adding together said blocks of magnitude- squared data to produce a peak, the 
location of said peak being determined using digital interpolation methods and 
corresponding to said pseudorange. 

36. The tracking process of claim 35 wherein said matched filtering 
operation comprises: 

performing a convolution of the compressed block's data against 
the pseudorandom sequence (PRS) of the GPS satellite being processed, said 
convolution being performed using said fast convolution algorithms to produce a 
product of the convolution. 

37. The tracking process of claim 36 wherein the fast convolution 
algorithm used in processing the buffered GPS signals is a Fast Fourier Transform 
(FFT) and the product of the convolution is produced by computing the forward 
transform of said compressed block by a prestored representation of the forward 
transform of the PRS to produce a first result and then performing an inverse 
transformation of said first result to recover said product. 

38. A computer readable medium containing a computer program 
having executable code for a GPS receiver, said computer program comprising: 

first instructions for receiving GPS signals from in view satellites, 
said GPS signals comprising pseudorandom (PN) codes; 

second instructions for digitizing said GPS signals at a 
predetermined rate to produce sampled GPS signals; 

third instructions for storing said sampled GPS signals in a 

memory; and 

fourth instructions for processing said sampled GPS signals by 
performing fast convolution operations on said sampled GPS signals, said fourth 
instructions comprising a matched filtering operation to determine the relative 
timing between said PN codes and locally generated PN reference signals. 

39. A method as in claim 12 wherein said GPS signals were sampled at 
a rate of a multiple of 1.024 MHz to provide said sampled GPS signals. 
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40. A GPS receiver as in claim 1 wherein said predetermined rate is a 
multiple of 1.024 MHz. 

41. A GPS receiver as in claim 1 wherein said DSP also performs 
preprocessing and postprocessing operations. 

42. A GPS receiver as in claim 41 wherein said preprocessing 
operation occurs before said fast convolution and said postprocessing operation 
occurs after said fast convolution. 

43. A GPS receiver as in claim 42 wherein said preprocessing 
operation comprises correcting for Doppler shift of signals from said in view 
satellite. 

44. A GPS receiver as in claim 42 wherein said preprocessing 
operation comprises adding together portions of said sampled IF GPS signals to 
provide compressed samples and wherein said fast convolution comprises a 
convolution of said compressed samples. 

45. A GPS receiver as in claim 44 wherein said fast convolution 
produces a plurality of results and wherein said postprocessing operation 
comprises adding together said plurality of results. 

46. A GPS receiver as in claim 45 wherein said plurality of results 
comprises a plurality of squares of magnitudes. 

47. A process utilizing global positioning system (GPS) satellites for 
determining the position of a remote sensor, the process comprising the steps of: 

receiving and storing GPS signals at said remote sensor from a 
plurality of in view GPS satellites. 

computing in the sensor, pseudoranges utilizing said GPS signals, 
said computing comprising digital signal processing using fast convolution 
techniques on stored GPS signal; 

receiving a transmission of satellite data information comprising 
data representative of ephemeris for a plurality of satellites; 

computing position information in the sensor by using said satellite 
data information and said pseudoranges. 

48. A process as in claim 47 wherein said transmission comes from a 
basestation. 

49. A process as in claim 47 wherein said transmission comprises 
transmissions from said plurality of satellites. 
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50. A process as in claim 47 wherein said position information is 
transmitted to a basestation. 

51. A process as in claim 48 further comprising receiving a precision 
carrier frequency signal from said basestation; automatically locking to said 
precision canier frequency signal from said basestation; and calibrating a local 
oscillator in said remote sensor with said precision carrier frequency signal. 

52. A process as in claim 49 wherein said remote sensor comprises a 
GPS receiver which receives said transmissions comprising data representative of 
ephemeris for a plurality of satellites. 

53. A process as in claim 47 wherein the step of computing the 
pseudoranges further comprises: 

storing the received GPS signals in a memory; 

processing the stored GPS signals for one or more of the in view 
GPS satellites in a digital signal processor by: 

breaking the stored data into a series of contiguous blocks whose 
durations are equal to a multiple of the frame period of the pseudorandom (PN) 
codes contained within the GPS signals; 

for each block, creating a compressed block of data with length 
equal to the duration of a pseudorandom code period by coherently adding together 
successive subblocks of data, said subblocks having duration equal to one PN 
frame; 

for each compressed block, performing a matched filtering 
operation to determine the relative timing between the received PN code contained 
within the block of data and a locally generated PN reference signal, said matched 
filtering operation utilizing said fast convolution techniques; 

determining said pseudorange by performing a magnitude-squared 
operation on the products created from said matched filtering operation and 
combining said magnitude-squared data for all blocks into a single block of data by 
adding together said blocks of magnitude-squared data to produce a peak, the 
location of said peak being determined using digital interpolation methods and 
corresponding to said pseudorange. 

54. The process of claim 53 wherein said matched filtering operation 
comprises: 
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performing a convolution of the compressed block's data against 
the pseudorandom sequence (PRS) of the GPS satellite being processed, said 
convolution being performed using said fast convolution algorithms to produce a 
product of the convolution. 

55. The process of claim 54 wherein the fast convolution algorithm 
used in processing the buffered GPS signals is a Fast Fourier Transform (FFT) 
and the product of the convolution is produced by computing the forward 
transform of said compressed block by a prestored representation of the forward 
transform of the PRS to produce a first result and then performing an inverse 
transformation of said first result to recover said product. 

56. A process as in claim 47 wherein said step of computing further 
comprises performing a preprocessing operation before said fast convolution 
techniques and performing a post processing operation after said fast convolution 
techniques. 

57. A process as in claim 56 wherein said fast convolution techniques 
comprise a matched filtering operation and wherein said GPS signals are stored in 
a series of contiguous blocks in a memory and wherein said preprocessing 
comprises, for each block, creating a compressed block of data by adding together 
successive subblocks of data, and wherein said post processing comprises adding 
together a representation of the products created from said matched filtering 
operation. 

58. A process as in claim 34 further comprising receiving a precision 
carrier frequency signal from said basestation; automatically locking to said 
precision canier frequency signal from said basestation; and calibrating a local 
oscillator in said remote sensor with said precision carrier frequency signal. 

59. A process as in claim 34 wherein said step of computing further 
comprises performing a preprocessing operation before said fast convolution 
techniques and performing a post processing operation after said fast convolution 
techniques. 

60. A process as in claim 59 wherein said fast convolution techniques 
comprise a matched filtering operation and wherein said GPS signals are stored in 
a series of contiguous blocks in a memory and wherein said preprocessing 
comprises, for each block, creating a compressed block of data by adding together 
successive subblocks of data, and wherein said post processing comprises adding 
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together a representation of the products created from said matched filtering 
operation. 

61. A GPS receiver as in claim 7 wherein said satellite data information 
comprises an identification of a plurality of satellites in view of said GPS receiver 
and a corresponding plurality of Doppler information for each satellite of said 
plurality of satellites in view of said GPS receiver. 

62. A GPS receiver as in claim 7 further comprising a power 
management circuit coupled to said downconverter and to said digitizer, wherein 
after said IF GPS signals are stored in said memory, said power management 
circuit reduces power consumed by said downconveter and said digitizer. 

63. A method of power management for a GPS receiver, and method 
comprising the steps of: 

receiving at said GPS receiver, GPS signals from in view satellites; 

buffering said GPS signals; and 

reducing power consumed by said GPS receiver. 

64. A method as in claim 63 further comprising: 

processing said GPS signals in a processing system to provide 
processed GPS signals. 

65 . A method as in claim 64 wherein said processed GPS signals 
comprise a pseudorange information. 

66. A method as in claim 65 further comprising: 

receiving at said GPS receiver a Doppler information of a satellite 
in view of said GPS receiver. 

67 . A method as in claim 66 further comprising: 
transmitting said pseudorange information. 

68. A method as in claim 64 wherein said step of reducing power 
consumed by GPS receiver occurs after said buffering of said GPS signals. 

69. A method as in claim 64 wherein said processed GPS signals are 
used to provide a latitude and a longitude of said GPS receiver. 

70. A method as in claim 64 wherein said.processing comprises a fast 
convolution operation on said GPS signals which are buffered in a memory. 

71. A method as in claim 70 wherein a digital signal processor (DSP) 
performs said processing. 
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72. A method as in claim 70 wherein the amount of said GPS signals 
which are buffered in said memory may be varied to trade off sensitivity with 
power reduction. 

73. A method as in claim 72 wherein less of said GPS signals may be 
buffered in order to save more power. 

74. A method as in claim 66 wherein said processing comprises a fast 
convolution operation on said GPS signals which are buffered in a memory. 

75. A method as in claim 74 wherein a digital signal processor (DSP) 
performs said processing. 

76. A method as in claim 74 wherein the amount of said GPS signals 
which are buffered in said memory may be varied to trade off sensitivity with 
power reduction. 

77. A method as in claim 76 wherein less of said GPS signals may be 
buffered in order to save more power. 

78. A method as in claim 66, wherein a communication receiver 
receives said Doppler information and said method further comprising reducing 
power consumed by said communication receiver for a period of time after 
receiving said Doppler information. 

79. A method as in claim 78 wherein said period of time is 
predetermined. 

80. A method as in claim 65 further comprising reducing power 
consumed by a transmitter which transmits said pseudorange information after 
transmitting said pseudorange information. 

81. A method as in claim 64 wherein a communication receiver receives 
satellite data information, and wherein said method further comprises: reducing 
power consumed by said communication receiver for a period of time after 
receiving said satellite data information. 

82. A method as in claim 81 wherein said satellite data information 
comprises data representative of ephemeris for a satellite. 

83 . A method as in claim 82 wherein said period of time is 
predetermined. 

84. A method as in claim 82 further comprising: 
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transmitting a position information by a communication transmitter 
in said GPS receiver, wherein power consumed by said communication transmitter 
is reduced after transmitting said position information. 

85. A method as in claim 84 wherein said position information 
comprises a latitude and longitude. 

86. A method as in claim 64 further comprising transmitting a position 
information by a communication transmitter in said GPS receiver, and reducing 
power consumed by said communication transmitter after transmitting said 
position information. 

87. A method as in claim 64 further comprising reducing power 
consumed by said processing system after providing said processed GPS signals. 

88. A method of managing power consumed by a system having a 
GPS receiver, said method comprising: 

receiving at said GPS receiver, GPS signals from in view satellites; 

receiving at a communication receiver coupled to said GPS receiver 
a signal comprising data representative of satellite data information; 

transmitting from a communication transmitter information 
representative of position information; 

processing said GPS signals in a processing system coupled to said 

GPS receiver, 

reducing power consumed by a component selected from the group 
consisting of said GPS receiver, said communication receiver and said 
communication transmitter. 

89. A method as in claim 88 wherein power consumed by said 
processing system is reduced. 

90. A method as in claim 88 wherein said step of reducing power 
comprises placing said component in one of an off state or a low power state. 

91 . A method as in claim 90 wherein said low power state comprises 
clocking said component at a lower clock frequency than when power to said 
component is not reduced. 

92. A method as in claim 88 further comprising: 

buffering said GPS signals in a memory and wherein said satellite 
data information comprises a Doppler information of a satellite in view of said 
GPS receiver. 
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93. A method as in claim 88 wherein after said step of reducing power, 
said component is in a reduced power state, said method further comprising 
returning said component to a normal power state upon receiving a signal through 
said communication receiver. 

94. A method as in claim 93 wherein said component is said 
communication receiver, said step of reducing power reduces power for said 
communication receiver for a period of time. 

95. A method as in claim 94 wherein said period of time is 
predetermined. 

96. A system in a GPS mobile unit having a reduced power state, said 
system comprising: 

a GPS receiver for receiving GPS signals; 

a communication receiver coupled to said GPS receiver, said 
communication receiver for receiving a signal comprising satellite data information; 

a processing system coupled to said GPS receiver, said processing 
system for processing said GPS signals; 

a communication transmitter coupled to said processing system, 
said communication transmitter for transmitting information representative of 
position information for said GPS mobile unit; and 

a power management circuit coupled to said communication 
receiver, said power management circuit for reducing power consumed by a 
component selected from the group consisting of said GPS receiver, said 
communication receiver, said processing system, and said communication 
transmitter, wherein said component is placed in a reduced power state. 

97. A GPS mobile unit having a reduced power state, said GPS mobile 
unit comprising: 

a receiver for receiving GPS signals from in view satellites; 

a memory coupled to said receiver to store data representing said 

GPS signals; 

a processor coupled to said memory, said processor processing 
said GPS signals to provide processed GPS signals; 

a power management circuit coupled to said processor, said power 
management circuit reducing power consumed by said GPS mobile unit. 
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98. A GPS mobile unit as in claim 97 further comprising a 
communication receiver and a communication transmitter being coupled to said 
power management circuit. 

99. A GPS mobile unit as in claim 97 wherein said power management 
circuit reduces power consumed by said processor. 

100. A system as in claim 96 wherein after said component is in said 
reduced power state, said power management circuit returns said component to a 
normal power state after receiving a signal from said communication receiver. 

101. A GPS mobile unit as in claim 96 wherein after said GPS mobile 
unit is in a reduced power state, said power management circuit returns said GPS 
mobile unit to an increased power consumption state after receiving a signal from 
said communication receiver. 

102. A GPS mobile unit as in claim 97 further comprising a battery and 
a solar cell and a power regulator coupled to said battery and to said solar cell and 
to said power management circuit, said solar cell for charging said battery. 

103. A GPS mobile unit as in claim 102 further comprising a second 
battery coupled to said power regulator, wherein when said solar cell charges said 
battery, said second battery provides power to said GPS mobile unit. 

1 04. A GPS mobile unit as in claim 97 further comprising: 

a first controlled power interconnect coupled to said receiver and to 
said power management circuit; 

a second controlled power interconnect coupled to said memory 
and to said power management circuit, wherein said power management circuit 
reduces power by controlling the power provided to said receiver over said first 
controlled power interconnect and by controlling the power provided to said 
memory over said second controlled power interconnect. 

105. A GPS mobile unit as in claim 97 wherein said power management 
circuit comprises a microprocessor and a plurality of power switches. 

106. A GPS mobile unit as in claim 97 wherein said power management 
circuit comprises power control logic in a digital signal processing component and 
wherein said processor comprises said digital signal processing component. 

107. A GPS mobile unit as in claim 106 wherein said power 
management circuit further comprises a plurality of power switches coupled to said 
power control logic. 
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108. A system as in claim 96 further comprising a battery and a solar cell 
and a power regulator coupled to said battery and to said solar cell and to said 
power management circuit. 

1 09. A GPS mobile unit as in claim 97 further comprising a 
communication receiver which receives a satellite data information which 
comprises a Doppler information of a satellite in view of said GPS receiver. 

110. A GPS mobile unit as in claim 97 wherein after said GPS mobile 
unit is in a reduced power state, said power management circuit returns said GPS 
mobile unit to an increased power consumption state. 

111. A GPS mobile unit as in claim 97 further comprising a 
communication receiver which receives a satellite data information which 
comprises data representative of ephemeris for a satellite. 

112. A system as in claim 96 further comprising: 

a first controlled power interconnect coupled to said GPS receiver 
and to said power management circuit; 

wherein said power management circuit reduces power by 
controlling the power provided to said GPS receiver over said first controlled 
power interconnect 

113. A system as in claim 96 wherein said power management circuit 
comprises a microprocessor and a plurality of power switches. 

114. A system as in claim 96 wherein said power management circuit 
comprises power control logic in a digital signal processing component and 
wherein said processor comprises said digital signal processing component. 

115. A system as in claim 114 wherein said power management circuit 
further comprises a plurality of power switches coupled to said power control 
logic. 

116. A system as in claim 96 wherein said satellite data information 
comprises a Doppler information of a satellite in view of said GPS receiver. 

1 17. * A system as in claim 96 further comprising: 

a buffer coupled to said processing system, said buffer for storing 
said GPS signals. 

118. A system as in claim 9697 wherein said satellite data information 
comprises data representative of ephemeris for a satellite. 
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ll 9. A GPS mobile unit as in claim 97 wherein said processor processes 
said GPS signals by performing a fast convolution operation on a representation of 
said GPS signals. 

1 20. A GPS mobile unit as in claim 1 19 wherein after storing said GPS 
signals in said memory, power consumed by said receiver is reduced. 

121. A GPS mobile unit as in claim 1 19 wherein a preprocessing 
operation on said GPS signals is performed before performing said fast 

convolution operation. 

1 22. A GPS mobile unit as in claim 1 21 wherein a postprocessing 
operation on the results of said fast convolution operation is performed after 
performing said fast convolution operation. 

1 23. A system in a GPS mobile unit having a reduced power state, said 

system comprising: 

a GPS receiver for receiving GPS signals; 

a communication receiver coupled to said GPS receiver, said 
communication receiver for receiving a signal comprising satellite data informauon; 

a processing system coupled to said GPS receiver, said processing 

system for processing said GPS signals; 

a communication transmitter coupled to said processing system, 
said communication transmitter for transmitting information representative of 
position information for said GPS mobile unit; and 

a power management circuit coupled to a component, said power 
management circuit for reducing power consumed by said component selected 
from the group consisting of said GPS receiver, said communicauon receiver, said 
processing system, and said communication transmitter, wherein said component 
is placed in a reduced power state by said power management circuit. 

124 A system as in claim 123 wherein said power management circuit 
comprises a timer, said timer providing a periodic timing signal to said power 
management circuit to provide a full power state to said component. 

125 A system as in claim 61 wherein said component is in said reduced 
power state for a predetermined period of time and is then returned to a full power 
state. 

1 26. A system as in claim 1 23 further comprising: 
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a buffer coupled to said processing system, said buffer storing said 

GPS signals. 

1 27. A system as in claim 126 further comprising: 

a first controlled power interconnect coupled to said component and 
to said power management circuit, said power management circuit reducing power 
by controlling the power provided to said component over said first controlled 
power interconnect. 

1 28. A system in claim 126 wherein said power management circuit 
comprises a microprocessor and a plurality of power switches. 

1 29. A system as in claim 126 wherein said power management circuit 
comprises power control logic in a digital signal processing component and 
wherein said processing system comprises said digital signal processing 
components. 

1 30. A system as in claim 129 wherein said power management circuit 
further comprises a plurality of power switches coupled to said power control 
logic. 

131. A system as in claim 126 wherein said satellite data information 
comprises a Doppler information of a satellite in view of said GPS receiver. 

1 32. A system as in claim 126 wherein said satellite data information 
comprises data representative of ephemeris for a satellite. 

1 33. A system as in claim 126 wherein said processing system 
processes said GPS signals by performing a fast convolution operation on a 
representation of said GPS signals. 

1 34. A system as in claim 133 wherein after storing said GPS signals in 
said buffer, power consumed by said communication receiver is reduced. 

1 35. A system as in claim 133 wherein a preprocessing operation on 
said GPS signals is performed before performing said fast convolution operation. 

1 36. A system as in claim 1 35 wherein a postprocessing operation on 
the results of said fast convolution operation is performed after performing said 
fast convolution operation. 

1 37. A GPS mobile unit as in claim 97 further comprising a local 
oscillator coupled to said receiver, said local oscillator providing a first reference 
signal and a communication receiver coupled to said local oscillator, said 
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communication receiver providing a precision carrier frequency signal for 
calibrating said local oscillator which is used to acquire said GPS signals. 

138. A system as in claim 131 further comprising a local oscillator 
coupled to said communication receiver, said local oscillator providing a first 
reference signal, said communication receiver coupled to said local oscillator, said 
communication receiver providing a precision carrier frequency signal for 
calibrating said local oscillator which is used to acquire said GPS signals. 

1 39. In a method for determining the position of a remote unit, a process 
comprising: 

receiving, at said remote unit, a Doppler information of a satellite in 
view of said remote unit; and 

computing, in said remote unit, a position information for said 
satellite by using said Doppler information. 

1 40. A process as in claim 1 39, further comprising: 
transmitting said Doppler information of said satellite from a 

basestation to said remote unit. 

141. A process as in claim 140 wherein said Doppler information is 
obtained from a GPS receiver at said basestation. 

142. A process as in claim 140 wherein said position information 
comprises pseudoranges to a plurality of satellites in view of said remote unit, 
including said satellite. 

1 43. A process as in claim 140 wherein said position information 
comprises a latitude and longitude which indicates the position of said remote unit. 

1 44. A process as in claim 142 further comprising: 
transmitting said pseudoranges from said remote unit to said 

basestation, and wherein said basestation computes a latitude and longitude which 
indicates the position of said rembte unit. 

1 45. A process as in claim 142 further comprising transmitting satellite 
data information of said satellite to said remote unit, said satellite data information 
comprising data representative of ephemeris for said satellite. 

1 46. A process as in claim 143 further comprising transmitting satellite 
data information of said satellite to said remote unit, said satellite data information 
comprising data representative of ephemeris for said satellite. 
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1 47. A mobile unit which uses data representative of GPS signals to 
provide the position of said mobile unit, said mobile unit comprising; 

a receiver in said mobile unit, said receiver operable for coupling 
through a communications link to receive a Doppler information of a satellite in 
view of said mobile unit; 

a processing unit in said mobile unit, said processing unit coupled 
to said receiver to receive said Doppler information and compute a position 
information for said satellite by using said Doppler information. 

148. A method of using a basestation for providing a communications 
link to a mobile GPS unit, said method comprising: 

determining a Doppler information of a satellite in view of said 
mobile GPS unit; 

transmitting said Doppler information of said satellite in view to 
said mobile GPS unit. 

149. A method as in claim 148, wherein said Doppler information 
represents the Doppler shift of GPS signals from said satellite to said basestation. 

1 50. A method as in claim 149 wherein said Doppler information 
approximately represents the Doppler shift of GPS signals from said satellite to 
said mobile GPS unit. 

151. A method as in claim 148 wherein said Doppler information is 
obtained from a GPS receiver at said basestation and wherein said Doppler 
information represents the Doppler shift of GPS signals from said satellite to said 
basestation. 

152. A method as in claim 151 wherein said Doppler information 
approximately represents the Doppler shift of GPS signals from said satellite 
mobile GPS unit. 

153. A method as in claim 152 further comprising: 

receiving a position information from said mobile GPS unit, said 
position information being received at said basestation such that said basestation 
obtains a latitude and longitude which indicates the position of said mobile GPS 
unit. 

154. A method as in claim 153 wherein said position information 
comprises pseudoranges to a plurality of satellites in view of said mobile GPS 
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unit, including said satellite, and wherein said basestation computes said latitude 
and longitude from said pseudoranges. 

1 55. A method as in claim 153 wherein said position information 
comprises said latitude and longitude. 

1 56. A method as in claim 148 further comprising: 

transmitting satellite data information of said satellite to said mobile 
GPS unit, said satellite data information comprising data representative of 
ephemeris for said satellite. 

1 57. A process as in claim 139 wherein a processing unit uses said 
Doppler information to compensate for a Doppler shift of GPS signals from said 
satellite. 

1 58. A mobile unit as in claim 147 wherein said processing unit uses 
said Doppler information to compensate for a Doppler shift of GPS signals from 
said satellite. 

159. A mobile unit as in claim 158 wherein said communication link 
comprises a radio frequency communication medium. 

1 60. A mobile unit as in claim 158 further comprising; 

a transmitter coupled to said processing unit, said transmitter for 
transmitting said position information. 

161. A mobile unit as in claim 160 wherein said position information 
comprises pseudorange to a plurality of satellites in view of said mobile unit. 

1 62. A mobile unit as in claim 160 wherein said position information 
comprises a latitude and longitude which indicates the position of said mobile unit. 

163. A mobile unit as in claim 158 wherein said processing unit 
comprises a digital signal processing integrated circuit (DSP) and wherein said 
DSP processes said GPS signals and said Doppler information using a fast 
convolution algorithm. 

1 64. A mobile unit as in claim 163 further comprising: 

a transmitter coupled to said processing unit, said transmitter for 
transmitting said position information. 

165. A mobile unit as in claim 147 wherein said receiver is operable to 
receive satellite data information of said satellite from a source other than said 
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satellite, wherein said satellite data information comprises data representative of 
ephemeris for said satellite. 

1 66. A basestation for providing a communication link to a mobile GPS 
unit, said basestation comprising: 

a source of a Doppler information of a satellite in view of said 
mobile GPS unit; 

a transmitter coupled to said source of said Doppler information, 
said transmitter for transmitting through said communications link said Doppler 
information to said mobile GPS unit. 

167. A basestation as in claim 166 wherein said source of said Doppler 
information is a storage unit coupled to said basestation, said storage unit storing 
precalculated approximate Doppler information for said satellite. 

1 68. A basestation as in claim 166 further comprising: 

a receiver for receiving a position information from said mobile 

GPS unit; 

a processor coupled to said receiver. 

169. A basestation as in claim 166, wherein said Doppler information 
represents the Doppler shift of GPS signals from said satellite to said basestation. 

170. A basestation as in claim 169, wherein said Doppler information 
approximately represents the Doppler shift of GPS signals from said satellite to 
said mobile GPS unit. 

171. A basestation as in claim 166, wherein said Doppler information is 
obtained from said source which comprises a GPS receiver at said basestation and 
wherein said Doppler information represents the Doppler shift of GPS signals 
from said satellite to said basestation. 

1 72. A basestation as in claim 171 wherein said Doppler information 
approximately represents the Doppler shift of GPS signals from said satellite to 
said mobile GPS unit. 

173. A basestation as in claim 168, wherein said position information is 
received at said basestation such that said basestation obtains a latitude and 
longitude which indicates the position of said mobile GPS unit. 

174. A basestation as in claim 173, wherein said position information 
comprises pseudoranges to a plurality of satellites in view of said mobile GPS 
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unit, including said satellite, and wherein said processor of said basestation 
computes said latitude and longitude from said pseudorange. 

175. A basestation as in claim 173, wherein said position information 
comprises said latitude and longitude. 

1 76. A basestation as in claim 166 wherein said transmitter is also for 
transmitting satellite data information of said satellite to said mobile GPS unit, said 
satellite data information comprising data representative of ephemeris for said 
satellite. 

177. A basestation as in claim 169 wherein said basestation and said 
mobile GPS unit are within approximately 150 kilometers of each other. 

178. A method of calibrating a local oscillator in a mobile GPS receiver, 
said method comprising: 

receiving a precision carrier frequency signal from a source 
providing said precision carrier frequency signal; 

automatically locking to said precision carrier frequency signal and 
providing a reference signal ; 

calibrating said local oscillator with said reference signal, said local 
oscillator being used to acquire GPS signals. 

1 79. A method as in claim 178 wherein said receiving step comprises 
extracting said precision earner frequency signal from a data signal containing 
satellite data information communicated over a communication link. 

1 80. A method as in claim 179 wherein said satellite data information 
comprises a Doppler information of a satellite in view of said mobile GPS 
receiver. 

181. A method as in claim 179 wherein said satellite data information 
comprises data representative of ephemeris for a satellite. 

1 82. A method as in claim 179 wherein said communication link is 
selected from the group consisting of a two-way pager link or a cellular telephone 
link or personal communication system or specialized mobile radio or a wireless 
packet data system. 

183. A method as in claim 179 wherein said communication link is a 
radio frequency communication medium. 
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1 84. A method as in claim 178 wherein said automatic frequency control 
logic comprises one of a phase lock loop or a frequency lock loop or a block phase 
estimator. 

1 85. A method as in claim 1 84 wherein said reference signal provides a 
reference frequency which is compared to a frequency provided by said local 
oscillator to calibrate said local oscillator. 

1 86. A mobile GPS receiver comprising: 

a first antenna for receiving GPS signals; 

a downcon verier coupled to said antenna, said antenna providing 
said GPS signals to said downconverter; 

a local oscillator coupled to said downconverter, said local 
oscillator providing a first reference signal to said downconverter to convert said 
GPS signals from a first frequency to a second frequency; 

a second antenna for receiving a precision canier frequency signal 
from a source providing said precision carrier frequency signal; 

an automatic frequency control (AFC) circuit coupled to said 
second antenna, said AFC circuit providing a second reference signal to said local 
oscillator to calibrate said first reference signal of said local oscillator, wherein said 
local oscillator is used to acquire said GPS signals. 

187. A mobile GPS receiver as in claim 1 86 further comprising a 
comparator coupled to said AFC circuit and to said local oscillator, said 
comparator comparing said first reference signal and said second reference signal 
to adjust the frequency of said first reference signal from said local oscillator. 

1 88. A mobile GPS receiver as in claim 1 87 wherein said AFC circuit 
comprises a phase lock loop coupled to a receiver which is coupled to said second 
antenna. 

1 89. A mobile GPS receiver as in claim 1 86 further comprising a 
receiver coupled to said second antenna, said receiver for receiving said precision 
carrier frequency signal from said second antenna, wherein said receiver receives 
said precision carrier frequency signal with a data signal containing satellite data 
information communicated through said second antenna. 

1 90. A mobile GPS receiver as in claim 1 89 wherein said satellite data 
information comprises a Doppler information of a satellite in view of said mobile 
GPS receiver. 
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191. A mobile GPS receiver as in claim 190 wherein said satellite data 
information comprises an identification of a plurality of satellites in view of said 
mobile GPS receiver and a corresponding plurality of Doppler information for 
each satellite of said plurality of satellites in view of said mobile GPS receiver. 

1 92. A mobile GPS receiver as in claim 1 89 wherein said satellite data 
information comprises data representative of ephemeris for a satellite. 

1 93. A method of using a basestation to calibrate a local oscillator in a 
mobile GPS receiver, said method comprising: 

producing a first reference signal having a precision frequency; 

modulating said first reference signal with a data signal to provide a 
precision carrier frequency signal; 

transmitting said precision carrier frequency signal to said mobile 
GPS receiver, said precision carrier frequency signal being used to calibrate a local 
oscillator in said mobile GPS receiver, said local oscillator being used to acquire 
GPS signals. 

1 94. A method as in claim 193 wherein said data signal contains satellite 
data information which comprises a Doppler information of a satellite in view of 
said mobile GPS receiver. 

1 95. A method as in claim 193 wherein said data signal contains satellite 
data information which comprises data representative of ephemeris for a satellite. 

1 96. A method for determining the position of a remote unit, comprising 
the steps of : 

transmitting GPS satellite information, including Doppler, to a 
remote unit from a basestation via a data link; 

receiving at said remote unit said satellite information and GPS 
signals from in view satellites; 

computing, in the remote unit, pseudoranges to said in view 

satellites; 

transmitting said pseudoranges to said basestation from said remote 
unit via said data link; and 

computing, at said basestation, the position of said remote unit, 
using said pseudoranges. 
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1 97. A basestation for providing a calibration signal for use in a mobile 
GPS receiver to calibrate a local oscillator in said mobile GPS receiver, said 
basestation comprising 

a first source for a first reference signal have a precision frequency; 

a modulator coupled to said first source and to a second source of 
satellite data information said modulator providing a precision carrier frequency 
signal; 

a transmitter coupled to said modulator, said transminer for 
transmitting said precision carrier frequency signal to said mobile GPS receiver, 
said precision frequency signal being used to calibrate said local oscillator, said 
local oscillator being used to acquire said GPS signals. 

1 98. A basestation as in claim 197 wherein said satellite data information 
comprises a Doppler information of a satellite in view of said mobile GPS 
receiver. 

199. A basestation as in claim 197 wherein said satellite data information 
comprises data representative of ephemeris for a satellite in view of said mobile 
GPS receiver. 

200. A basestation as in claim 197 further comprising a processor 
coupled to said transmitter, said processor instructing said transmitter to transmit 
to said mobile GPS receiver. 

201 . A basestation as in claim 200 wherein said processor determines a 
plurality, of satellites in view of said mobile GPS receiver and obtains said satellite 
data information for each satellite of said plurality of satellites, and wherein said 
processor instructs said transmitter to transmit to said mobile GPS receiver an 
identification of said plurality of satellites and said satellite data information. 

202. A basestation as in claim 201 wherein said satellite data information 
comprises Doppler information for said plurality of satellites. 

203. A basestation as in claim 201 wherein said satellite data information 
comprises data representative of ephemeris for said plurality of satellites. 

204. A process as in claim 143 wherein said position information further 
comprises altitude of said remote unit. 
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FIGURE 5B 
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